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In bacteria such as Escherichia coli, DNA is compacted into a nucleoid near the cell center, 
while ribosomes—molecular complexes that translate messenger RNAs (mRNAs) into proteins—are 
mainly localized to the poles. We study the impact of this spatial organization using a minimal 
reaction-diffusion model for the cellular transcriptional-translational machinery. While genome-wide 
mRNA-nucleoid segregation still lacks experimental validation, our model predicts that ~ 90% of 
mRNAs are segregated to the poles. In addition, our analysis reveals a “circulation” of ribosomes 
driven by the flux of mRNAs, from synthesis in the nucleoid to degradation at the poles. We show 
that our results are robust with respect to multiple, biologically relevant factors, such as mRNA 
degradation by RNase enzymes, different phases of the cell division cycle and growth rates, and the 
existence of non-specific, transient interactions between ribosomes and mRNAs. Finally, we confirm 
that the observed nucleoid size stems from a balance between the forces that the chromosome 
and mRNAs exert on each other. This suggests a potential global feedback circuit in which gene 
expression feeds back on itself via nucleoid compaction. 


I. INTRODUCTION 

The cytoplasm of many bacterial cells exhibits a strik¬ 
ing spatial organization: rather than filling the entire cell 
volume, the DNA forms a condensed structure called a 
“nucleoid” that is generally localized near midcell (Fig. 1) 
[1, 2]. Moreover, ribosomes—large molecular complexes 
that translate messenger RNAs (mRNAs) into proteins— 
are observed to be anti-localized from the nucleoid [2]. 
These observations raise two natural questions: (1) What 
physical processes are responsible for this subcellular or¬ 
ganization? (2) How does this internal structure influence 
the basic processes of mRNA transcription and protein 
translation in the cell? 

In the model bacterium Escherichia coli {E. coli), ~ 
1.5 mm of supercoiled DNA are compacted into a ~ 1 /rm^ 
nucleoid volume [3], thus forming a dense DNA mesh with 
average pore diameter ~ 50 nm. As a result, free ribo¬ 
somes, with diameter ^ 20 nm, can readily diffuse into 
the nucleoid [4], while polysomes, molecular complexes 
composed of mRNAs with multiple bound ribosomes and 
having an effective diameter > 50 nm, anti-localize from 
the nucleoid due to excluded-volume effects. In vivo mea¬ 
surements of mRNA mobility suggest a typical diffusion 
coefficient of D ~ 0.05 /xm^/s, implying that mRNAs 
formed in the nucleoid by transcription from DNA can 
diffuse out of the nucleoid to the ribosome-rich regions 
in a few seconds—a time significantly shorter than the 
typical mRNA lifetime of ~ 5 min [2, 5]. These observa¬ 
tions suggest that most mRNAs formed in the nucleoid 
diffuse out of the nucleoid to the ribosome-rich regions 
where ribosomes and mRNAs are colocalized, and where 
the bulk of translation occurs. 

In recent years, the idea that mRNA localization may 
play a functional role in bacteria [6, 7] has inspired a vari¬ 
ety of measurements of mRNA localization. These stud¬ 
ies have provided evidence for multiple, mRNA-speciflc 
localization patterns, such as localization within the cell 


cytoplasm [8], to the cell membrane or to the cell poles 
[9], and at the nascent septum separating daughter cells 
[7]. However, such experiments have proven to be chal¬ 
lenging [10, 11] and limited to specific mRNAs [12, 13[. 
As a result, genome-wide ribosome-mRNA colocalization 
still lacks experimental validation. In this study, we em¬ 
ploy known reaction-diffusion properties of mRNAs, ribo¬ 
somes, and the nucleoid to predict the physical origin and 
the extent of overall, genome-wide mRNA localization. 

Our approach also describes non-specific, transient in¬ 
teractions between ribosomes and mRNAs. In this re¬ 
gard, recent studies in E. coli [2[ and Caulohacter cres- 
centus [14[ observed an increase of the diffusion coefficient 
of non-translating ribosomes under depletion of the pool 
of mRNAs, which has been interpreted as evidence for 
non-specific, transient bindings between ribosomes and 
mRNA molecules in [14[. On the other hand, another 
study in E. coli reported that the diffusion coefficient 
of non-translating ribosomes is not affected by mRNA 
depletion [4[. Here, we show that if transient ribosome- 
mRNA interactions are significantly faster than other rel¬ 
evant time scales, the reaction-diffusion equations can be 
substantially simplified by treating such interactions as a 
local Poisson process. Finally, we consider the opposing 
forces exerted by the compressed nucleoid and by polarly 
localized mRNAs and confirm that the observed nucleoid 
size results from the balance of these dominant forces, 
and we present a simple analytical formula for the degree 
of nucleoid compaction under different physiological con¬ 
ditions, e.g. for different amounts of DNA and mRNA in 
the cell. 


II. RESULTS 

We describe the coupled dynamics of ribosomes and 
mRNAs in an E. coli cell using a minimal, ID reaction- 
diffusion model. We introduce a coordinate x running 
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FIG. 1: Schematic of the spatial organization of transcrip¬ 
tion and translation in E. coli. mRNAs (black solid curves) 
are transcribed in the nucleoid—the condensed DNA chromo¬ 
some at the cell center (gray solid curve). Ribosomes (red 
circles) bind to mRNAs forming polysomes, i.e. mRNAs with 
multiple bound ribosomes. Polysomes diffuse preferentially 
out of the nucleoid due to excluded-volume effects. Eventu¬ 
ally the mRNA molecules are degraded (black dashed lines). 
In our ID model, the coordinate x runs along the long axis 
of the cell and, assuming symmetry, we model only the right 
half of the cell {0 < x < £). 


along the long axis of the cell and, given the approx¬ 
imate left-right symmetry of a typical E. coli cell, we 
consider only the right half, from a; = 0 at the cell center 
to a; = I! at the right cell pole (Fig. 1). In vitro mea¬ 
surements of the assembly dynamics of the translation- 
initiation complex suggest that the in vivo binding rate of 
308 ribosomal subunits to mRNAs is significantly larger 
than the unbinding rate, thus implying that the majority 
of mRNAs have a 308 subunit bound at the translation- 
initiation site [15]. If so, translation is largely governed 
by the dynamics of 508 ribosomal subunits, and therefore 
we initially consider only the 508 units, which we refer to 
simply as “ribosomes”. We further assume that ribosomes 
may undergo transient, non-specific binding to mRNAs; 
extended versions of our model including the two ribo¬ 
somal subunits and disallowing non-specific ribosome- 
mRNA interactions are presented later, and they con¬ 
firm qualitatively the results obtained with the simple 
model discussed here—see 8upporting Information, sec¬ 
tions 810 and 814, for details. The ID concentration of 
free (F) ribosomes, Cf.{x), denotes the number of F ri¬ 
bosomes per unit length in an infinitesimal slice of the 
cell perpendicular to the x axis. 8imilarly, Pm,nix) is the 
ID concentration of mRNAs with m transiently bound 
(B) [14] ribosomes and n translating (T) ribosomes. As 
shown in section 81, the average number of ribosomes 
per mRNA, m -|- n ~ 12, is well below the maximum to¬ 
tal number of ribosomes, m + n ^ 100, that could be 
linearly packed onto a typical mRNA. Thus, we consider 
only mRNA species with m < rumax) n < Umux, where 
Wniax and riniax are some maximal numbers of allowed ri¬ 
bosomes per mRNA, chosen large enough to account for 
all typical mRNA species present in the cell. Importantly, 
this choice reduces the number of mRNA species present 
in our model, thus making it computationally tractable. 

The resulting reaction-diffusion equation for the F- 
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where only mRNA species with allowed values of 0 < m < 
m-tnax and 0 < n < n„iax are considered. The first term 
on the right-hand side (RH8) represents diffusion includ¬ 
ing excluded-volume effects due to the condensed DNA 
(Fig. 1). In this term, Dp is the diffusion coefficient for F 
ribosomes, which incorporates crowding effects due to ri¬ 
bosomes, mRNAs and other macromolecules, while Vp{x) 
is the fractional volume available to an F ribosome within 
the DNA mesh at position x (Fig. 81), see sections 82 
and 84 for details. The first term in the second line repre¬ 
sents F ribosomes binding to all possible mRNA species 
with m B ribosomes and n T ribosomes, and thus becom¬ 
ing B ribosomes. This term is proportional to the F —>■ B 
transition rate and to the total density of mRNA. 
In principle, should decrease with the ribosome oc¬ 
cupancy number m + m however, here m -1- n is much 
smaller than the maximum packing density, thus this ef¬ 
fect is small, see above. The next term in the second line 
describes a B ribosome unbinding from an mRNA, where 
denotes the unbinding rate and the multiplicity fac¬ 
tor m accounts for multiple B ribosomes on the mRNA. 
8imilarly, the third line represents transitions between F 
and T ribosomes, where is assumed independent of 
m, n (see above). 8ince measurements suggest that the 
lifetime of the B state is significantly shorter than that of 
the T state [14], here the F —B —T transition is incor¬ 
porated into the F —^ T transition, with an effective rate 
Finally, the last line represents B and T ribosomes 
being freed from mRNA molecules as these are degraded 
at rate (3 [16]. 

The ID reaction-diffusion equations for the mRNA 
densities are 
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Here D is the average mRNA diffusion coefficient in the 
cytoplasm, and Vm+nix) is the fractional available volume 
within the nucleoid for an mRNA with m + n attached 
ribosomes (Fig. 81), see section 84 for details. The third 
and fourth lines represent binding and unbinding of B 
















3 



FIG. 2: DNA fluorescence along the long cell axis for 3 fj,m- 
long E. coll cells grown in glucose minimal media. Cells 
were stained with SYTOX Orange and imaged at exponen¬ 
tial phase. Fluorescence for a few representative cells (gray), 
and resulting average fluorescence over 35 cells with standard 
error of the mean (red), both symmetrized and normalized to 
unit area. Left inset: nucleoid vs. cell length, where the gray 
area includes cells with length with 5% of 3 pm. Right inset: 
ribosomal protein S2-YFP (green) and nucleoid (red). Scale 
bar: 1 pm. 


and T ribosomes from mRNAs of species to, n, while 
the fifth line represents B ribosomes binding to mRNAs 
of species to — 1, n or unbinding from mRNAs of species 
TO-l-1, n. Similarly, the terms in the sixth line represent T 
ribosomes binding and unbinding from mRNAs of species 
TO,n — 1 and m,n + 1. Finally, the last line represents 
transcription of initially ribosome-free mRNAs according 
to the nucleoid profile a{x), and mRNA degradation at 
rate (3. 

As transcriptional and translational time scales 
(< 1 min) are fast compared to cell doubling times 
(>20min), we focus on steady-state conditions. At 
steady state, in Eq. (1) we enforce a constraint 
on the total number of ribosomes, 2 dx\c^{x) -f 
+ 'n)Pm,u{x)] = ^tot, and we 

set a no-flux boundary condition at the cell pole, 
[dcj,{x)/dxv^(x) — Cj,{x)dv-p{x)/dx\\^^i = 0, see section 
S2. Similarly, in Eq. (2) we impose no-flux bound¬ 
ary conditions at the cell pole and at the cell center, 
[dpm,n{x)/dxVm+n{x)-Pm,n{x)dVm+n{.x)ldx%=o^l = 0 , 
the latter reflecting the left-right symmetry of the cell. 
According to this symmetry, the flux of F ribosomes at 
midcell must also vanish, and this follows directly from 
the boundary conditions above—see section S5. 

We fix the model parameters from experimental data 
as follows. We considered E. coli cells in glucose minimal 
media with a ^ 1/h growth rate, see section S3, and fo¬ 
cused on the mid-phase of the division cycle by selecting 


cells with length within 5% of a typical, medium length 
of 2^ = 3pm, compare left inset in Fig. 2. We then 
rescaled the resulting DNA fluorescence profiles to a cell 
length of 2 £, and we estimated the nucleoid profile along 
the long cell axis by averaging over multiple cells, see 
main panel in Fig. 2 and section S4. The F ribosome dif¬ 
fusion coefficient was taken to be Dp = 0.4pm^/s, while 
the diffusion coefficient of mRNAs was set at the aver¬ 
age diffusion coefficient of polysomes, D = 0.05pm^/s 
[2, 4]. These diffusion coefficients were obtained from 
measurements of mean square displacements of free and 
mRNA-bound ribosomal subunits in living E. coli cells 
[2], and thus include the effects of cytoplasmic crowding. 
The F-ribosome available volume Vp(x) was estimated 
from the concentration profiles of free ribosomal subunits 
[4], and similarly for the mRNA available-volume profiles 
Vm+n(x) (section S4). The mRNA degradation rate was 
taken to be /3 = 3 x 10“^/s, corresponding to a mean 
mRNA lifetime of ~ 5 min [5]. The total mRNA pro¬ 
duction rate Otot = 2 dxa{x) was obtained from the 
total number of mRNAs per cell via the relation atot//? = 
A^mRNA = 5 X 10^ [17], while the profile of mRNA produc¬ 
tion a{x) was chosen to be proportional to the DNA den¬ 
sity ^{x), see section S4 for details. The average time for 
a ribosome to complete protein translation is estimated 
to be Tt = 40 s [18], so we took the T —>■ F transition 
rate to be = 1 /tt = 2.5 x 10“^/s. To set the other 
transition rates we used the observation that ~ 80% of 
ribosomes are T ribosomes [2], with non-T ribosomes es¬ 
timated to spend ~ 90% of their time as B ribosomes and 
~ 10% of their time as F ribosomes, where the 90%: 10% 
division of non-T ribosomes between B and F ribosomes 
is inferred from ribosome diffusion rates in Caulohacter 
crescentus [14]. Global equilibrium between T and F ri¬ 
bosomes at steady state requires the F —>■ T transition 
rate to satisfy Nppk'^^ = where Np and 

are the total number of F and T ribosomes, respectively, 
and p ~ 1.7 x 10^/pm is the average total mRNA ax¬ 
ial density [17]. Also, in the equilibrium condition above 
we neglected the rate PNp at which mRNA-bound T ri¬ 
bosomes are freed by mRNAs that are being degraded, 
because this rate is much smaller than Npk'^g. According 
to the estimate above for the average number of F, T, 
and B ribosomes, we have N^/Np = 80%/(20% x 10%), 
which yields k'^^ = 6 x lO^'^pm/s. The total number 
of ribosome in the cell was taken to be Wot = 6 x 10'^ 
[2[. As for the the binding-unbinding rates k^^, fc^g of 
the F -fA B transition, the equilibrium condition reads 
Np p k^^ = Np k^g where, since the F o B transition 
occurs on timescales not longer than Is [14], we ne¬ 
glected the niRNA-decay term jdNp <C Np k^g. To¬ 
gether with the estimation above for the average num¬ 
ber of F, T, and B ribosomes Np/Np = 2%/18%, this 
equilibrium condition provides an estimate for the ratio 
kpp/k^g ~ 5.4 X 10“^ pm, but it does not specify the indi¬ 
vidual values of k^^, k^g. However, because the timescale 
of the F o B transition is significantly shorter than other 
relevant timescales [14], these processes can be treated in 
the limit where they are at rapid equilibrium. In this 
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FIG. 3: Steady-state mRNA and polysome distributions. To¬ 
tal mRNA density ptot(*) (red) and density pm,n.{x) of mR- 
NAs with m transiently bound (B) ribosomes and n trans¬ 
lating (T) ribosomes for 0 < m < mmax = 8 and 0 < n < 
nmax = 24 (gray). The density po,o(x) of ribosome-free mR- 
NAs (green) and the density P 8 , 24 (x) of mRNAs with the 
largest number of T and B ribosomes considered (blue) are 
also shown. The profiles Pm,n.(x), pop(x), and p 8 . 24 (a:) are 
normalized to unit area. Inset: distribution of mRNA species, 
shown as a heat map of the number Nm,n of mRNAs with m 
B ribosomes and n T ribosomes in the right half of the cell. 
The maximal number of T ribosomes per mRNA used in our 
model, Umax = 24, is indicated. 


limit, the problem can be significantly simplified, and 
the set of 1 -b (mmax + l)(''T-max + 1) reaction-diffusion 
Eqs. (1), (2) reduces to a set of nmax + 2 equations which 
completely characterize the solution Cp(x) and Pm,n{x) 
for any m and 0 < n < rimax- Importantly, these rapid- 
equilibrium equations depend on only through 

their ratio—see section S6 for details. 

We numerically solved Eqs. (1), (2) at steady state 
in the rapid-equilibrium limit for B ribosomes, by fixing 
the maximal number of allowed T ribosomes per mRNA 
at rimax = 24, see section S6 for details. The resulting 
mRNA profiles and ribosome concentrations are shown 
in Figs. 3 and 4, respectively. As shown in Fig. 3, the 
total mRNA profile ptot{x) = Em=o 
markedly localized outside the nucleoid region—most of 
the mRNAs are segregated from the nucleoid because of 
excluded volume. Indeed, the density profiles Pm,n{x) 
show that the more ribosomes an mRNA has bound, 
the more segregated the mRNA is from the nucleoid, see 
Fig. S2 for details. Since mRNAs are generated at the 
nucleoid, the strong segregation of mRNAs away from 
the nucleoid at steady state must result from the ma¬ 
jority of mRNAs being loaded with ribosomes so that 
excluded volume biases their diffusion away from the 
nucleoid. This conclusion is confirmed in the inset of 
Fig. 3, which shows a heat map of the total number 
Nm,n = Jq dxpm,n{x) of niRNAs of species m,n. Most 
mRNAs are loaded with ~ 10 T ribosomes and ^ 2 B ribo¬ 
somes. These average loading numbers are in agreement 
with the experimental estimates above for the ribosome 
numbers [2, 17]: The number of T ribosomes per mRNA 
is A^t/A^hirna ~ (80% Wot)/iVmRNA ~ 10, and a similar 
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FIG. 4: Steady-state ribosome concentrations and ribosome 
fluxes. Top: Goncentrations Ct{x), Cb{x), and Cf(®) of trans¬ 
lating (T), transiently bound (B), and free (F) ribosomes in 
the right half of the cell (compare Fig. 1). Bottom: Fluxes 
of T, B, and F ribosomes along the cell’s long axis depicted 
in the top panel. The arrow length is proportional to local 
ribosome flux, and the arrows in the legends correspond to a 
flux of 20/s. 


estimate yields ^2 B ribosomes per mRNA. In addition, 
the inset shows that the chosen value rimax = 24 is large 
enough to encompass all typical mRNA species that are 
present. 

Since each ribosome has a linear size of ~ 20 nm [19], 
the effective size of an mRNA molecule with ^ 10 bound 
ribosomes is significantly larger than the pore size of the 
DNA mesh in the nucleoid, which we estimate to be ~ 
50 nm, see section S4 for details. Thus, the majority of 
mRNAs experience strong excluded-volume effects which 
push them out of the nucleoid region. 

Since mRNAs are created by transcription in the nu¬ 
cleoid but end up segregated away from the nucleoid, 
there must be a flux of mRNAs toward the cell poles. 
Given that new mRNAs are rapidly loaded with T ri¬ 
bosomes at a rate k^^Np/{2i) = 6 x 10“^/rm/s x 
1200/3p,m ~ 0.24/s, implying full occupation by ^ 10 
T ribosomes in ~ 3 s, the poleward flux of mRNAs car¬ 
ries with it a poleward flux of ribosomes. Since ribo¬ 
somes are conserved in our model, reflecting the long 
half-life of ribosomal components ]20], there must be a 
compensating flux of F ribosomes from the poles toward 
the nucleoid. In Fig. 4 we show the steady-state con¬ 
centration of F ribosomes CF(a;), the concentrations of 
T and B ribosomes, c^{x) = ^ Pm,n{x), 

Cb{x) = the ffux Jp = 

—Dp[dcp(x)/dxvp(x) — Cp(x) dvp(x)/dx] of F ribosomes, 
and the fluxes and Jb of T and B ribosomes, compare 
Eqs. (S45), (S46). As expected from the observed segre¬ 
gation of mRNAs, the T and B ribosomes are markedly 
excluded from the nucleoid, and there is a net poleward 
ffux of T and B ribosomes. Notably, the effect of ex¬ 
cluded volume in the nucleoid is so strong that mRNAs 
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and their associated ribosomes flow from a low- to a high- 
concentration region. By contrast, F ribosomes are small 
enough to freely penetrate the nucleoid, and a flux of F ri¬ 
bosomes is driven by the gradient of these ribosomes from 
the poles to the nucleoid. Overall, these results illustrate 
and quantify a “circular” flux for the ensemble of T, B, 
and F ribosomes, compare Fig. 1: First, multiple F ribo¬ 
somes bind to mRNAs made in the cell nucleoid. Each 
mRNA is thus loaded with ~ 10 T ribosomes and ~ 2 B 
ribosomes to become a polysome. Second, the effects of 
excluded volume in the nucleoid result in a net flow of 
these polysomes to the cell poles. Once polysomes reach 
the poles, they ultimately decay and free their ribosomes. 
This “pumping” of T and B ribosomes from the nucleoid 
to the poles results in an excess of F ribosomes at the 
poles, and thus in a diffusive return flux of F ribosomes 
to the nucleoid. 

The existence of a steady ribosome circulation implies 
that there must be an external source of energy driving 
these circular fluxes. There are two possible candidates 
within our model: Process (A) is the non-equilibrium 
creation and degradation of mRNAs, and Process (B) is 
mRNA and F-ribosome binding in the nucleoid and sub¬ 
sequent expulsion from the nucleoid by excluded-volume 
effects. Process (A) should be strictly dependent on new 
mRNA production, whereas Process (B) should persist 
even in the limit where the mRNA production and degra¬ 
dation rates are both low, with the total number of mR¬ 
NAs fixed and equal to Nm^NA- Therefore, we varied 
the mRNA production and degradation rates together, 
keeping the total mRNA number constant: The circula¬ 
tion vanished as the mRNA rates slowed (Fig. S3), thus 
identifying Process (A), the flux of new mRNAs from nu¬ 
cleoid to pole, as the driver of ribosome circulation. This 
conclusion is confirmed by an analytical estimate for the 
poleward flux of T and B ribosomes, which is shown to be 
proportional to the mRNA production rate, see section 
S8 for details. 

Before discussing other implications of our results, it 
is worth considering that mRNA transcription takes a 
finite amount of time. For the average mRNA length 
of ~ 3 X 10^ nt discussed in section SI and an average 
mRNA elongation speed of ~ 50nt/s [21], we obtain 
a typical transcription time of ~ 1 min, during which 
nascent mRNAs are bound to DNA while being elon¬ 
gated. We therefore extended our model to include these 
nucleoid-bound mRNAs, whose axial densities we denote 
by ni^)- these mRNAs turn into free mRNAs at a rate 
7 = 1/(1 min), see section S9 for details. Besides confirm¬ 
ing the picture in the model with one mRNA species, this 
extended model gives novel insights into the mechanism 
of co-transcriptional translation, namely the observation 
that ribosomes translate nascent mRNAs while these are 
being transcribed in the nucleoid [22]. While it has been 
previously hypothesized that most of the ribosomes in the 
dense nucleoid region translate co-transcriptionally [2], 
our analysis shows that only ^ 34 % of these ribosomes 
carry out co-transcriptional translation, while a compa¬ 
rable fraction of ~ 37 % translates post-transcriptionally. 
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FIG. 5: Steady-state ribosomal-subunit concentrations and 
fluxes for the model including SOS and 50S ribosomal subunits. 
Top: Concentrations C 7 os(a:), Csosb]*), C 3 osF(a;), c'^osix) of 
70S, transiently bound (B) SOS, free (F) SOS subunits, and of 
SOS subunits bound to the translation-initiation site, respec¬ 
tively. We also show the concentrations Csosb]®) and Csosf]®) 
of B and F SOS subunits, respectively. Bottom: Fluxes 
Jjosix), Jiossix), J 3 osf(x) of 70S, B and F SOS subunits, and 
flux J 30 S (a;) of SOS subunits bound to the translation-initiation 
site. The fluxes J^ossix), J 5 osf(x) of B and F SOS subunits 
are also shown. Fluxes are represented along the cell’s long 
axis depicted in the top panel, the arrow length is propor¬ 
tional to local ribosome flux, and the arrows in the legends 
correspond to a flux of SO/s. 


mostly on polysomes loaded with a relatively small num¬ 
ber of ribosomes. The extended model also allows us 
to address the effects of co-transcriptional translation 
on the protein-synthesis rate: introducing the efficiency 
£ = Jo dxp„,^n(x)/JVtot, i.e. the av¬ 

erage number of proteins translated per unit time per 
ribosome, we And that co-transcriptional translation im¬ 
plies a ~ 3% increase in ribosome efficiency, under the 
assumption of B ribosome binding to all transcripts, see 
section S9 for details. 

We then extended our model to take account of both 
ribosomal subunits. During translation initiation, a 30S 
subunit binds to the mRNA initiation site first, and then 
a 50S subunit binds to the 30S subunit forming a trans¬ 
lating 30S-50S (70S) pair [23]. To model this process 
and the spatiotemporal dynamics of the two ribosomal 
subunits, we introduce mRNA species with I B 30S sub¬ 
units, m B SOS subunits, n 70S pairs, and no 30S sub¬ 
unit at the initiation site, and denote their density by 
Pi,m,n(x,t). Similarly, we denote by p] ^ „(a^)i) the den¬ 
sity of mRNAs with I B 30S subunits, m B SOS subunits, 
n 70S pairs, and a 30S subunit at the initiation site. As- 
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suming rapid equilibrium for B 308 and 508 subunits, 
we solved the reaction-diffusion equations for the mRNA 
densities pi^m,n{x), concentrations 

CsosFix,t), C 5 osFix,t) of free 308 and 508 subunits—see 
section 810 for details. The results shown in Figs. 5, 818, 
and 87 confirm the picture obtained with the simpler one- 
subunit model. In particular, the mRNAs, both with and 
without a 308 subunit on the translation-initiation site, 
are strongly segregated from the nucleoid, and the larger 
I, m, n, the stronger the segregation. In Fig. 5 we show 
the concentrations and fluxes of the ribosomal subunits: 
the 708 subunits, the B 308 and 508 subunits, and the 
308 subunits bound to translation-initiation sites are all 
strongly segregated from the nucleoid, compare Fig. 4, 
while F 308 and 508 subunits are only slightly excluded 
from the nucleoid. Note that the concentrations of F and 
B 308 subunits are lower than the corresponding concen¬ 
trations of F and B 508 subunits because for our choice 
of parameters the majority of non-translating 308 sub¬ 
units are bound to translation-initiation sites, thus reduc¬ 
ing the number of F and B 308 subunits. The fluxes of 
708, B 308 and 508 subunits, and of 308 subunits bound 
to translation-initiation sites are directed toward the cell 
poles, while the compensating fluxes of F 308 and 508 
subunits are directed toward the cell center. 

The two-subunit model was then extended to include 
additional, biologically relevant features: In section 811 
we incorporated in the model the mechanism of mRNA 
degradation by RNase enzymes [16, 24], see Figs. 89- 
812, and in section 812 we extended the model to cells in 
the late phase of their division cycle, see Figs. 814-816. 
We then extended the analysis above for glucose minimal 
media to different growth rates: We imaged E. coli cells 
in glycerol minimal and defined rich media, with growth 
rates ~ 0.5/h and ~ 2/h, respectively, extracted the nu¬ 
cleoid profiles, and present the resulting model predic¬ 
tions in section 813, compare Figs. 817-820. Finally, as 
the existence of B ribosomes is an open question [2, 4, 14], 
in section 814 we considered a version of the model with 
two ribosomal subunits and no B subunits. Overall, these 
results confirm all the qualitative features of the simpler, 
two-subunit model discussed above. 

8o far, our analysis has shown that the excluded vol¬ 
ume due to DNA localization at midcell segregates the 
majority of mRNAs to the cell poles. In what follows, 
we will extend this analysis and show that the converse 
is also true, i.e. mRNA segregation to the poles causes 
nucleoid compaction at midcell. 8pecifically, in section 
815 we show that mRNAs, like particles in a gas, ex¬ 
ert an entropic force (pressure) on the nucleoid directed 
towards the cell center, and that this force can be com¬ 
puted directly from the reaction-diffusion Eqs. (1) and 
(2). On the other hand, the natural tendency of the 
compressed DNA polymer to increase its configurational 
entropy results in a an effective “spring” force pushing 
outward on the mRNAs. Exploiting the condition that 
these two forces must balance at mechanical equilibrium, 
we self-consistently determined the mRNA and F ribo¬ 
some profiles, as well as the DNA density profile ifiix), 


and the results are shown in Figs. 823 and 824. In par¬ 
ticular, denoting by 2xo the width of the nucleoid region 
centered at midcell, the resulting value of 2xo « 1.43/rm 
matches well the experimentally observed nucleoid size in 
Fig. 2. Finally, a simplified version of our analysis pro¬ 
vides a straightforward prediction for the nucleoid width 
under different physiological conditions: We assume that 
mRNAs are uniformly distributed outside the nucleoid, 
which they cannot penetrate, that the nucleoid is con¬ 
fined in a region of width 2xo centered at midcell, and 
we neglect the small force exerted by F ribosomes on the 
nucleoid (see section 815 for details). As a result, the nu¬ 
cleoid size xq can be determined by solving the following 
force-balance equation: 

iVmRNA TT^ 

2(£-xo) 6 (27ri?2)2/3a,5/3 ’ 

where ^ = 200 nm is twice the estimated persistence 
length of a segment of supercoiled DNA, N = 7.5 x 10^ 
is the total number of such segments in the chromosome 
[25], and R = 0.5 p,m is the radius of a circular cell slice. 
In Fig. 825 we show the predicted nucleoid size for dif¬ 
ferent values of the total number of mRNAs—the larger 
the number of mRNAs, the more the nucleoid shrinks to¬ 
wards midcell due to the entropic force exerted by the 
polysomes. 


III. DISCUSSION 

The study of intracellular mRNA localization has at¬ 
tracted growing interest in recent years [10-13]. In eu¬ 
karyotic systems, mRNA localization is a well-established 
mechanism for achieving a variety of functions, such as 
the targeted expression of proteins to specific regions of 
the cell, the control of intracellular signaling, or the par¬ 
tition of mRNAs into daughter cells for cell-fate differen¬ 
tiation [26, 27]. 8ome functional, mRNA-specific local¬ 
ization patterns has also been observed in bacteria [6]: 
Two examples are the membrane localization of mRNAs 
which code for proteins transporting lactose into the cell 
[9] and mRNA localization at the cell poles, which has 
been shown to play a functional role in controlling sugar- 
utilization genes [9, 28]. 

In this study, we analyze the extent of bacterial 
genome-wide mRNA localization by means of a minimal 
reaction-diffusion model for the transcriptional and trans¬ 
lational machinery in E. coli. The experimental observa¬ 
tion that ribosomes in E. coli are markedly localized to 
the cell poles, and thus segregated from the nucleoid lo¬ 
cated at the cell center ]2], suggests that the majority 
of mRNAs are also likely to be segregated from the nu¬ 
cleoid. While experiments on mRNA localization in E. 
coli have proven to be challenging [10, 11], and are limited 
to specific mRNAs [12, 13], our model makes a novel pre¬ 
diction for strong, genome-wide mRNA localization away 
from the nucleoid, indicating that ~ 90% of mRNAs are 
typically located outside the nucleoid, and demonstrating 
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that the total mRNA profile resembles that of translat¬ 
ing ribosomes (compare Figs. 3, 4). A specific predic¬ 
tion of our model is that mRNA segregation is due to 
excluded-volume effects resulting from the condensed nu¬ 
cleoid DNA. Overall, this result provides novel insights 
into the mechanisms of mRNA segregation: While other 
studies for specific mRNAs raised the possibility that 
mRNA segregation may be associated with dynamical re¬ 
arrangements of the nucleoid [29], our analysis indicates 
that genome-wide mRNA segregation can arise entirely 
from excluded-volume effects. Also, our result that seg¬ 
regation primarily affects mRNAs loaded with multiple 
ribosomes is in line with the recent experimental obser¬ 
vation that mRNAs with multiple bound fluorescent pro¬ 
teins localize to the cell poles in live E. coli cells [10]. 

Our model also reveals a “circulation” of ribosomes 
within the cell driven by the flux of newly synthesized 
mRNAs from the nucleoid to the poles: mRNA-bound 
ribosomes flow from the nucleoid to the cell poles, where 
they unbind from mRNAs and then diffuse back to the 
nucleoid to bind newly synthesized mRNA molecules. 
Using our model, we also analyze the extent of co- 
transcriptional translation, namely the observation that 
ribosomes translate mRNAs that are being transcribed 
in the nucleoid [22]. Although it has been recently 
hypothesized that most of the ribosomes in the DNA- 
rich region translate co-transcriptionally [2], we find that 
only about a third of these ribosomes carry out co- 
transcriptional translation, whereas a slightly higher frac¬ 
tion translates post-transcriptionally in polysomes with a 
relatively small loading number. 

We incorporated in our model the mechanism of mRNA 
degradation by RNase enzymes, different phases of the 
cell division cycle, different growth rates, and the effect 
of non-specific, transient interactions between ribosomes 
and mRNAs [2, 4, 14], showing that our results are stable 
with respect to such variety of conditions. 

Finally, we extended our analysis to study the conse¬ 
quences of mRNA localization on nucleoid compaction. 
Using our calculated ribosome and mRNA densities, we 
confirmed that mRNA segregation to the poles quantita¬ 
tively accounts for nucleoid compaction. Physically, the 
observed nucleoid size reflects the balance of two com¬ 
peting entropic forces—the compressive force that mR¬ 
NAs (polysomes) at the poles exert on the DNA, and 
the expansive force exerted by the DNA on these mR¬ 
NAs. Our detailed analysis supports a simplified ana¬ 
lytical formula (3) that predicts nucleoid size for differ¬ 
ent physiological conditions, compare Fig. S25. Biologi¬ 
cally, the compaction of the nucleoid by mRNAs creates a 
potential global feedback circuit: gene expression drives 
mRNA levels, which, by compacting the nucleoid, im¬ 
pact transcription-factor access and hence gene expres¬ 
sion [30]. 

To summarize, while localization of the transcriptional- 
translational machinery is a well-known, functional mech¬ 
anism in eukaryotes, the function of such localization in 
bacteria is not yet well established. In this regard, our 


model provides novel insights into the mechanisms gov¬ 
erning the spatial structure of transcription and transla¬ 
tion in bacteria, and can help guide the molecular manip¬ 
ulation of these functions, with potentially broad appli¬ 
cations in molecular synthetic biology and biotechnology. 
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SI Fraction of available ribosome binding sites 


In this section, we will estimate the fraction of available ribosome binding sites on an mRNA chain by making use 
of experimental data on mRNA synthesis and degradation in E. coli. To begin with, we estimate the average number 
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of nucleotides per mRNA, which we quantify as the ratio between the total number of nucleotides synthesized per unit 
time, and the total number of mRNAs synthesized per unit time. First, the number of nucleotides synthesized per unit 
time is given by the number of active RNA polymerases (RNAPs), i.e. ^ 1000 RNAP/cell, times the chain-elongation 
rate at which a single RNAP synthesizes an mRNA, i.e. ~ 50nt/s [1], resulting in a rate of ~ 5 x 10‘^nt/s. Second, 
the total number of mRNAs synthesized per unit time is given by the total number of mRNAs iVmRNA = 5 x 10^ 
[2] times the mRNA degradation rate /3 = 3 x 10“^/s [3], and it is given by ~ 15mRNA/s. The resulting average 
number of nucleotides per mRNA is ~ 3 x 10^ nt. Given that each ribosome covers ~ 30 nt [4], we obtain that there 
are on average ~ 3 x 10^ nt/(30nt) = 100 available ribosome binding spots per mRNA. Finally, we observe that there 
is a total number of iVtot ~ 6 x 10^ ribosomes per cell [5], i.e. there are on average A^tot/AmRNA ~ 12 ribosomes 
per mRNA, which is only ~ 12% of the number of available binding spots per mRNA, thus implying that the vast 
majority of the ribosome binding sites are free. 


S2 Diffusion equations in the presence of excluded volume 


In this section we derive the diffusion equation for a population of F ribosomes, taking into account the effects of 
excluded volume due to the nucleoid. The continuum limit of this equation can be derived starting from a discrete 
version [6]. We divide the interval 0 < x < £, from midcell to the cell pole, into bins of width Ax, where Xi denotes the 
position of bin i: Here, each bin represents a projection on the one-dimensional x axis of a three-dimensional section 
of the cell obtained by slicing the cell perpendicularly to its long axis. The master equation for the F ribosome density 
in bin i, c^{xi,t), can then be obtained by considering hopping of ribosomes between adjacent bins: 


dcp{xi,t) 

Ft 


dp {[c^{xi-i,t) + Cp{xi+i,t)]vt.{xi) - c^{xi,t)[v^{xi-i) + UF(xi+i)]} , 


(SI) 


where dp is the bare rate at which a ribosome hops from one site to another, and Up(x) is the fraction of available 
volume at position x, i.e. the probability that a ribosome finds enough free volume to hop into the DNA mesh at 
position X. The continuum limit of Eq. (SI) can be obtained by expanding Cp(x,t) and Up(x) around x = Xi, which 
yields 


9cp(x, t) 

Ft 


Df 


92 


Cp(x,t) (fvp{x) 


(S2) 


where Dp = dpAx^ is the bare diffusion coefficient. The ID F ribosome flux at position x is given by 


dp(x) — Tip 


dcp(x) , , , ,dxp(x) 

—-Up(x) - Cp(x)—-- 

ax ax 


(S3) 


Note that in the absence of reaction terms, the steady-state solution of Eq. (SI) is 

Cp(x) = C Up(x), 


(S4) 


where the constant C depends on the total number of F ribosomes. 


S3 E. coli nucleoid imaging analysis 

E. coli strain AFS55 [5] which contains ribosomal protein S2-YFP was inoculated overnight in the media including 
MOPS defined rich, glucose and glycerol minimal media (Teknova). Next day the overnight cultures were diluted by 
a hundred fold in the same fresh media and grew in 37 °C. 

The staining procedure is adopted from a previous publication [7]. Briefly, 100 /iL was taken when the OD reached 
0.3 and stained with 1/iL 50/iM SYTOX Orange (Molecular Probes) for 10 min at 37 °C in darkness. Cells were 
then washed twice with 1 mL fresh media and resuspended in 1 mL fresh media. 1 /iL of cells were placed on 1 % 
low-melting agar pad (Calbiochem) made from the same media and imaged with inverted NikonOOi epifluorescent 
microscope equipped with a 100 x 1.4 NA objective (Nikon) and Hamamazu Orca R2 CCD camera. NIS Elements 
software (Nikon) was used to automate image acquisition for phase contrast, YFP and mCherry fluorescent channels. 
Segmentation, quantification of fluorescence intensity, and cell-length measurements were further analyzed in MATLAB 
using customized programs. 
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S4 Available-volume profiles 


In this section we estimate the available volume for ribosomes, mRNAs, and polysomes inside the nucleoid. To 
estimate the DNA pore size inside the nucleoid, we recognize that the E. coli chromosome is characterized by a 
branched, plectonemic structure of supercoiled DNA [8]: We denote by L the total length of the plectonemic DNA 
polymer, and by V the volume wherein this polymer is confined. For simplicity, suppose that the DNA is arranged to 
occupy the edges of a three-dimensional cubic lattice of total volume V. The volume is divided into M cubic pores 
each with edge length a, so th at V = Ma^. In addition, the DNA length is related to the pore size by L = 3aM, thus 
providing the estimate a = ^/SV/L. For an estimated total plectoneme length L = 1.5 mm [8] and a nucleoid volume 
V = 1.2 /im^ [5], the estimated pore size is a ~ 50 nm. 

We now provide a rough, dimensional estimate of the available volume for an F ribosome inside the nucleoid. 
Continuing to use the simple cubic lattice model for the nucleoid, we estimate the excluded volume inside each cubic 
pore of DNA: as the ribosome approaches one of the lattice edges, the excluded volume is given by the edge length a 
times the area within the pore from which the center of mass of the ribosome is excluded when it approaches the edge, 
i.e. 1/4 of a circle with the ribosome’s radius r ~ lOnm [9]. Multiplying by the total number of edges (12), we obtain 
an estimate for the excluded volume, Vexci = 12a7rr^/4 = Sairr^. We then express the available-volume fraction in 
the nucleoid, which we will denote by in terms of the ratio between the excluded volume and the pore volume as 



The numerical correction factor n is introduced to improve our dimensional estimate, and is set by fitting our model 
to the experimentally observed concentration profile of freely diffusing ribosomes. If we imagine dividing the cell 
into two regions along the axial direction, the nucleoid region (in) and the polar region outside the nucleoid (out), 
then according to section S2 the concentration of freely diffusing ribosomes satisfies Since the 

concentration of freely diffusing ribosomes is, roughly speaking, about 10% larger outside the nucleoid than inside 
[10], we have ~ 90%, where = 1 because there is no DNA in the out region. As a result, we obtain 

K = 0.25 in Eq. (S5). 

The available-volume fraction for mRNAs loaded with m B ribosomes and n T ribosomes can be estimated along 
the same lines. Since the excluded volume for a single ribosome depends on its two-dimensional cross section 7rr^ 
rather than on its volume, we introduce an effective radius rm+n for an mRNA of species m, n, where Vm+n is chosen 
to reproduce the overall cross section of the mRNA along with its ribosomes: 

Trrl^+n=Trrl + {m + n)Trr^, (S6) 


and Tr = 20 nm is a typical mRNA radius of gyration [11]. The mRNA excluded volume in the nucleoid is then 
determined by replacing r with rm+n in Eq. (S5). 

Finally, the smooth excluded-volume profiles Vp(x), Vm+n{x) shown in Fig. SI were obtained from u™, Vm+n ns 
follows. First, we introduce a local density of DNA length 


ip{x) oc 


1 

1 -I- exp[((x/.^ - 1/2)] ’ 


(S7) 


with C = 20: this density profile is shown in Fig. SI, it was chosen to reproduce the experimental axial DNA 
fluorescence in Fig. 2, and normalized to the total DNA length per volume [5] by setting 


1 

£ 



dx (p{x) 


L 

V' 


(S8) 


Second, we replaced LjV in Eq. (S5) with ^{x), and we obtained 


and similarly 


Vp(x) = exp[—K7rr^(fi(x)], 

(S9) 

Vm+n(x) = exp[-K7rr^_^„(p(x)]. 

(SIO) 

Vr(x) = expl-Kirrlip^x)] 

(Sll) 


Denoting by 
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FIG. SI: Density of DNA length and available-volume profiles in ID model of the E. coli transcriptional-translational machinery. 
Density of DNA length ip{x) normalized to unit area, available-volume profile for free ribosomes Vp{x) = exp\—K-Kr^(p{x)], and 
available-volume profiles for polysomes with m transiently bound and n translating ribosomes Vm+nix) = exp[—K7rr^_|_„^(a;)] 
for right half of cell (compare Fig. 1), where the arrow indicates the direction of increasing m -I- n. Here the constant k = 0.25 
has been estimated by fitting the observed concentration profile of freely diffusing ribosomes [10], r = 10 nm is the ribosome 
radius [9], and Vm+n = \/’'b + (wi- + n)r'^, Tr = 20 nm are the polysome effective radius and the mRNA radius of gyration [11], 
respectively. 


the available volume for a ribosome-free mRNA, Eqs. (S6), (S9), (SIO) imply that the polysome available volume can 
be rewritten as 


Um+„(a;) = ■UR(a;)[uF(a;)]""+”. (S12) 

The factorized form of Eq. (S12) has a simple intuitive interpretation if we recall that available volumes Vm+nix) 
represents the probability that an mRNA of species m, n finds enough free volume to hop into the DNA mesh at 
position X, compare section S2. In this regard, Eq. (S12) treats the components of an mRNA of species m, n 
independently: the hopping probability Vm+n(x) is equal to the product of the hopping probability Vb,{x) for the bare 
mRNA and of the hopping probability Vp{x) for each of the m + n ribosomes carried by the mRNA. 


S5 Flux of free ribosomes at midcell 


In this section we show that at steady state the no-flux boundary conditions at the cell poles for F ribosomes and 
at the cell poles and midcell for niRNAs in Eqs. (1) and (2) imply that the flux of F ribosomes at midcell is also 
zero. This absence of F-ribosome flux at midcell is a necessary physical condition due to the assumptions of ribosome 
conservation and cell symmetry in our model. 

First, we use Eq. (S3) to rewrite the diffusive term in Eq. (1) at steady state as 


Dp 


d'^Cf.{x) 

dx'^ 


Vp{x) 


Cp{x) 


d^v^{x) 

dx^ 


dJj,{x) 

dx ’ 


(S13) 


we then integrate Eq. (1) at steady state with respect to x between a; = 0 and x = i, and we obtain 


i-e r 

Jp(0)-b / dx - fc®„Cp(a:)^^p™,„(a:)-f - fc;)’„Cp(a;)^^p^,„(a;)-bfc;Cff^^np,„,„(a:) (S14) 

m n m n m n m n 

+/3^ ^ (W -b n) Pm.nix) 


= 0 , 


where we used Eq. (S13) and the no-flux boundary condition at the cell pole, Jf{£) = 0, enforced in Eq. (1). We 


now consider Eq. (2) at steady state. We rewrite the diffusion term as 


dPm,n{^) 

dx 




dvn 




dx 
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multiply both sides by m, sum with respect to m and n, and integrate with respect to x. We obtain 


0 = 


/ dx'^ - k^^Cpix) ^ m Pm,n{x) - Pm,n{x) 

^ n ^ m—0 m—1 

^max '^max 1 '^max 

+ k^^C:p(x) mpm-i,nix) + k^fi ^ m{m + 1) Pm+i,nix) -/S rnpm,n{x) 


m—1 


m—0 


m—0 


'^max 1 


pt "‘'max "‘'max "‘max 

/ dx^ kl^C^{x) ^ Pm,n{x) - k^g ^ m Pm,n{x) - p ^ m Pm,n{x) 
^ n ^ m—0 m—1 m—0 


(S15) 


where we wrote explicitly the summation bound in terms of rumax) and we obtained the third line by changing 
summation variables in the first and second terms in the second line. Multiplying Eq. (2) by n, summing with respect 
to m and n, integrating with respect to x and proceeding along the same lines, we obtain 


^max 1 


pt “max “max “max 

/ dx^ kl^C^{x) ^ Pm,n{x) - ^ U p^^ri{x) - P ^ n Pm,n{x) 

m ^ n—0 n—1 n—0 


= 0. 


(S16) 


Combining Eqs. (S14), (S15), and (S16), we obtain Jf(0) = 0, i.e. there is no F-ribosome flux at midcell. 


S6 Rapid equilibrium for transiently bound ribosomes 


In this section, we consider the reaction-diffusion processes of ribosomes and mRNAs in the limit that transient 
ribosome binding and unbinding is significantly faster than all other relevant processes. As a result, the statistics of 
B ribosomes can be described entirely as a local Poisson process, without the need for a separate kinetic equation 
for each mRNA species with a different number of B ribosomes. This substantially simplifies the reaction-diffusion 
equations, with a major impact on computational tractability. 


Let us consider Eq. (2), the reaction-diffusion equation for mRNAs from the main text, at steady state: 

D 


d Pra^nipO / \ / \^ 

9 Pm.ny^) 9 

dx^ dx^ 


+ k^^C^{x)pjn-l,n{x) + k^ff{m+l) Pjn+l,n{x) + kl^C^{x)pm,n-l{x) + fcjff (n + 1) p^.n+l (a:) 


jd Pm,n{,x') — 0. (S17) 


In the rapid-equilibrium limit, where k^^, k^g are both large, we set 

=AC, fcoff = AfcBff, (S18) 

where A ^ 1 and k^g are of order unity. We now consider Eq. (2) to leading order in A, and we have 

^off l)Pm+l,n(a:) ^on^F(^)P7n,n (a") 0, 0 ^ ^ ^ ^max 1? 

which implies 

PmA^) = ^ j POA^) = Pm,nA)- (S20) 

Equation (S20) implies that the distribution of B ribosomes at position x is Poissonian: this is a direct consequence 
of the rapid-equilibrium limit, where at any position x binding and unbinding of B ribosomes dominates over all other 
processes, such as binding-unbinding of T ribosomes and diffusion. 

Now, we sum Eq. (S17) with respect to m = 0,... ,TOmax) and we obtain the following conservation equation for 
the density of mRNAs at position x 

'^max 

E 

m—0 

“h/pon^F (^)/^m,n—1 (^) ^off Pm,n-\-l{.^') r ^n,0 f^i^x'j ^ ^ ^ — 0- 

m—0 


D- 


dx 


dpm,n (^) 
dx 


'^m-\-n{^) Pm,ni^) 


dx 


- Knd 


’(a;)pm,„(a;) - k'^gup^^nix) 
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Equations (S20) and (S21) to leading order in A imply the following equation involving only po,?^ 

dJm,n (^) 


E 

m—O 


dx 




+ 6„fia{x) = 0, 


where we simplified the notation introducing the local flux of mRNAs of species m, n 


j / \ _ 7-) I dpm,nix) , , , .dVm+n{x) 


dx 


dx 


(S22) 


(S23) 


and defining 


— ^on^F ^off ^ (^) ~f ^on^F {^)P^Ti,n—1 (^) “f ^ofF ~f l)Pm,n+l(^) [d Pna^nis^^ • (S24) 


Using Eqs. (S12), (S20), and (S23), and setting 

pn{x) = Po,n{x), 

we compute explicitly the first term in Eq. (S22) for large m-aiax- 


(S25) 


E 


dx 






dx ^ V dx 

m—0 ^ 




m—0 

dVm+nix) 




m—0 


m—1 


J_ f konCi.{x) \"" fcg„ dc^jx) 7.1 / fcg„CF(a:) A”' dp„(a:) 
m!^ V k^ff ) k^ft dx ^ ml \ k^(r ) dx 


'■off 


'-off 




dx 

UR(a;)[r;F(x)]""+" 

\ 

m+n-l dVp{x) 




) bF(x)]^ 


dx 
dvp{x) 


kon / \ t \ ( dc^ix) 

r^Pn[x)Vji{x) —- V^{x} - Cj,{x)—— 

Y dx dx 


dp„{x) .dv^{x) / ^ / X 1 dv^{x) 

H-^-^'R(a;) - Pn{x)— - np„(x)v^{x}—^— - 

dx dx Vf.{x) dx 


(S26) 


The sums with respect to m in the second line of Eq. (S26) are all of the form J2m=o w™/to!, and we calculated 
them explicitly with a change of summation variables in some terms. The sum of ujm,n in Eq. (S22) can be computed 
along the same lines as in Eq. (S26). We substitute Eq. (S26) in Eq. (S22), and we obtain the final set of rimax + 1 
reaction-diffusion equations for p„: 


^Tx\ 


k^„c 


on».7F(x) 

UB 

'^off 


/ X \ r / xinl 7 X / X dC:p{x) dVj.{x)\ dp„{x) 

Uf(x) ) [uF(a:)] I -^Pn{x)vnix) “ CF(a:) j H- 


.dvj^{x) / X / X 1 dvp{x) 

-Pn{x)—, - np„{x)v^{x) 


dx 


v-p{x) dx 


+ exp [-A:;^„CF(a:)p„(x) - k’^gup^ix) 


+k'^^Cj.{x)pn-iix) + k'^ff{n+l)pn+iix)-I3pnix)]+6„fia(x) = 0. (S27) 


Next, we derive the rapid-equilibrium equation for the concentration of F ribosomes starting from the reaction- 
diffusion equation (1) in the main text. To achieve this, we introduce the subleading correction to Pm,n, setting 


Pm,n{x) — Pm,nix^ -f ^^Pm,n{x) • 

We then substitute Eq. (S28) into Eq. (S17), and we obtain: 

kosi'IXl + l)pm+l.„(x) - k^^Cj.{x)pm,n{x) = 


(S28) 


kosi'ixi + 1) 


Pm-t-l,n(^) “t” . ^Pm+l,n{x) 


A 


- konC^ix) 


Pm,nix^ T ^^Pm,nix') 


-\- 1) 7^(x) A^q^^Cf (x) y, (x) 


(S29) 
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where in the third line we observed that the terms of order A cancel out because of Eq. (S20). Importantly, Eq. (S29) 
implies that the term in the first line, which represents the current between mRNA species (m, n) and (to + 1, n), can 
be expressed in terms of the combination of subleading corrections to pm,n in the last line. 

Using Eqs. (S23), (S24), (S29) in Eq. (S17), we obtain the following set of equations 

[fc^fjTO Apm,n{x) - fcQ„Cp(a;)Ap„_i_„(a;)] - [k^s{m + l)Apm+i,n{x) - k^^Cp{x)Apm,n{x)] = 

_ dJm^n{x) (- 2 ,) _)_ SmfiSnfiCtix). (S30) 

ax 

Exploiting the recursive structure of Eqs. (S30) [12], we solve them iteratively for the quantities k^^m Apm,n{x) — 
k^^Cp{x)Apjn-i,n{x), and we obtain 


k^^mmaxApra^^^,n{x) k^^C-p{x)Apm^^^-l,n{x) — 
^off(n^max ~ ~ k^^Cp{x)Apm^^^-2,n{x) = 


dJm^^^,n{x) 

dx 

dJm^^^-l,n{x) 

dx 

dJm^^^,n{x) 

dx 


^rn^^^,n{x), 

^m„ax,n(^)j 


(S31) 


kof^ixi Apjji jii^x^ k^^^Cp(^x^Apiji—i ji(^x') — ^ ^ 


p—m 


dJp^nix) 

dx 


■UJp,nix) 


koS^PiAx) - fco„Cp(a;)Apo.n(a;) = ^ 

p^l 


dJp^nix) 

dx 


■UJp^nix) 


(S32) 


Now, let us compute the RHS of Eq. (1) to leading order in A. The second and third term on the RHS read 


^max 1 ^max 


‘'max “'max 


‘max '“max 


k^nC--eix) E E Pm, nix) + k^ff E E m Pm, nix) = E E 


m—Q n=0 


m—1 n—0 


n=0 m—1 

^max "^max ^^max 

= E E E 

n—O m—1 p—m 


^max "^max 

= E E 

n—O m—1 


m 


dx 

dJm,n ix) 
dx 


o,nix) 


■ ^m,nix) 


(S33) 


where in the first line we used Eq. (S29), in the second line we used Eqs. (S31), ..., (S32), and in the last line we 
rewrote the double sum over to, p as a single sum over to by using the identity 


^max '^max 


E E “p = E 

m—1 p—m m—1 


(S34) 


which is valid for any sequence am- The first term in the last line of Eq. (S33) can be computed for large TOmax along 



the lines of Eq. (S26): 


‘^max ^^max 

E E 

n—O m—1 


m 


^max OO 




n—0 m—1 


1 2 f konCF{x)\"' k^^dc^ix) 


-m 


l,B 

'‘■off 


J -Pnix) + —m 

ax ml 


1 (k^^c^{x)\^ dpn{x) 


pB 

'‘off 


dx 


_1 ( k^^c^jx) 

■ a-b 

'‘off 


Pn{x) 


dvFjjx) 

dx 


-"i EI -p 


n—0 


konCpix) 
pB 
'‘off 


vAx) bF(a;)]' 


Vj,{x)]^~^^ + Vnix){m + n)[vj.{x)] 


1 , fc?„CF(a;)„, fc^n dcj^ix) 

~r , X3 V-p^X) I , X3 T 


dJrn,n{,x) 

dx 

v^{x)[v,{x)r+'^ 

\ 

m+n-1 dv^jx) 


dx 


pB 

'‘off 


koS dx 


- Pn{x)v^{x)v^{x) 


k^^c^{x) dpn{x) .... k^^c^{x) .dv^{x) 

-Tb -h- V^{x)v^{x) --g- Pn{x)— - V^{x) 

k^^ dx k^f^ dx 


konCF^x) ( fconCF(a;) 


pB 

'‘off 


, fc^Cpfx) , ... .dvp{x) fc^Cpfx) , , , , 

1 H-- -Vf{x) Pn{x)Vt,{x)— -- -pn{x)vR(x)n 

koff J dx 


dVf.{x) 

dx 


The second term in the last line of Eq. (S33) reads 

^max ^^max 

mUJm,n{x) = -/? E E 1XiPm,n{x), 


"^max ^^max 


>. (S35) 


(S36) 


n=0 m—1 


n—0 m—1 


where we exchanged the sums with respect to n, m, and we observed that the only term in ujm,n that is nonzero when 
summed with respect to m and n is the term proportional to /3, see Eqs. (S24). The remaining terms in the RHS of 
Eq. (1) can be computed by using Eqs. (S20) and (S25), by which we calculate the exponential sums with respect to 
TO. Using Eqs. (1), (S33), (S35), and (S36), we obtain the final equation for the concentration of F ribosomes in the 
rapid-equilibrium limit 


E {-p 


n=0 


konCF^x) 


'‘off 


Vp{x) I [uF(a:)]’^ 




1 + 


d^cJx) , , , ,d'^v-p{x) 

Vp(Xj — Cp{x) 


'‘off 


dx'^ 

dCf{x) 

dx 


Vp{x) — Cp{x) 


dx^ 

dvF{x) 


y<^Pu{x)vAx) + 


fc®„CF(a;) / dpn{x) 


'‘off 


'‘off 


dx 


-Vr{x) - Pn{x) 


- exp 


konCFjx) 

kB 

'‘off 


dvFjjx) 

dx 

"^max 1 


,(a:) - 


kon^Fix) 


'‘off 


Pn{x)v^{x)n 


dx 

dvp{x) 

dx 


-kl^c^{x) Pn{x) + /Coff npn{x) +/3 npn{x) I = 0, (S37) 

n—0 n—1 n—0 / 

where the terms in the second and third line represent the diffusive flux of B ribosomes carried by diffusing mRNAs. 

Let us now discuss the boundary conditions for the rapid-equilibrium equations: The boundary conditions for Eq. 
(S27) are 




dx dx 

In Eq. (S37), the no-flux condition boundary condition at the right pole reads 

= 0 . 


= 0 . 


X—0,£ 


dcp{x) dv^ix) 

Vf.{x) — Cp{x) 


dx ' dx 

while the constraint on total ribosome number can be written as 


(S38) 


(S39) 


Ntot = 2 / dx 


= 2 dx 


Cf.{x) -f (to -I- n)pm,n{x) 


m—0 n—0 
uB 


konCF{x)\ ^ ( fc®„Cp(a;) 


»,(x) + exp X: 

V ''off / „=n 


.1=0 ^ 


Pn{x) 


(S40) 
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where in the second line we summed with respect to m by using Eqs. (S20) and (S25). 

Overall, the rimax + 2 reaction-diffusion equations (S27), (S37) and their boundary conditions (S38), (S39), (S40) 
completely characterize the solution Pnix), Cp{x). Using Eqs. (S20) and (S25), we obtain from Pnix) and Cp{x) the full 
set of mRNA concentrations Pm,n{x) for any m and 0 < n < rimax as well as all other physical quantities of interest, 
such as the total mRNA concentration 


Ptot(2;) — ^ ^ ' Pm,n(x) 


= exp 


—75- ) 

^off / „=o 


the concentrations of T and B ribosomes 


c.(x) = X! X! '^Pm,n{x) 


= exp 


k^^c^{x) 


^ npnix), 


i{x) = X! X! 


%ff V '‘-off / n=0 


the total mRNA flux 




= —Uexp 


[u,(x)]" ^p^{xW{x) (^^^V,{x) - cAx)^ 

^off / r)=n ^off N 


dpn{x) -.dv^ix) / N / N 1 dvT,{x) 

H- - - Vn{x)-Pn{x)— - npn{x)VR(x) yw —1- , 

ax ax Vf{x) ax 


and the fluxes of T and B ribosomes 




5; n[vAx)r\^Pu{x)v^ix) (^^vAx) - c^x)^^ 


= —Dexp 


dpn{x) .dv^ix) / N / X 1 dvp{x) 

H-1-Vb.(x) - Pn(x)—-np„(x)VR(x)—yw —1 - 

dx dx v-p(x) dx 


— ^ ^ ^ ^ dm,nix') 


-Dexp "£'i..wr [ (i + 

k^^ k^^c^ix) (dpnix) dvR{x)\ k^^c^{x) dv^{x) 

koS k°g \ dx dx J k°g dx 


The solution of the rapid-equilibrium equations allows for computing another important quantity, the ribosome 
efficiency, which we define as the average number of proteins translated per unit time per ribosome. We thus estimate 
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FIG. S2: Comparison of mRNA species inside versus outside the nucleoid. (A) Heat map of the number of mRNAs of species 
m, n inside the nucleoid „ Pm,n(x) in the right half of the cell as a function of the numbers m, n of transiently 

bound and translating ribosomes, respectively. (B) Same plot as in the left panel for the number of mRNAs outside the nucleoid 

= Si/2 Pra,ri{x). 


the total rate of protein production per cell in our model as the total rate at which T ribosomes fall of mRNAs after 
completing translation, i.e. 


S = 2kSs'^ X! / dxpm,n{x) 


m—0 n—1 


= 2 kSg n dx exp 


kon^Fix) 

k^ 


Pnix), 


and we define the ribosome efficiency as £ = S/Altot- 


(S47) 


S7 mRNA species inside versus outside the nucleoid 


In Fig. S2 we show the distribution of mRNA species in the nucleoid region and in the polar region as heat maps 
for the numbers N™ of mRNAs of species m, n for 0 < a: < ^|2 and Z/2 < x < I, respectively. 


S8 Estimate of ribosome flux 


In this section we present an analytical estimate of the flux of mRNA-bound ribosomes flowing from the nucleoid 
region toward one cell pole. The ribosome flux can be estimated by assuming that every newly synthesized mRNA 
migrates from the nucleoid to the pole, and carries an average load of ribosomes, yielding an estimated ribosome flux 


Oitot AA- + A^b 
2 AinRNA 


(S48) 


where atot = 2 dxa{x) ~ 17/s is the total mRNA transcription rate in the nucleoid region, and the second factor is 
the average number of ribosomes carried by each mRNA, with N^. = 4.8 x 10"^, = 1.1 x 10^, A/brna = 5 x 10^. 

Equation (S48) predicts a flow of J®+b ~ 100/s, in rough agreement with our numerical result Jt+b = Jt + As ~ 30/s 
for the flux of T and B ribosomes at the nucleoid edge shown in Fig. 4. The factor of ~ 3 discrepancy is not due 
to our estimate of the mRNA flux T^rna — Q^tot/2 — 8.5/s, which agrees reasonably well with the numerical result 
AnRNA — 7/s at the edge of the nucleoid. Rather, the difference is due to the fact that mRNAs with fewer ribosomes 
diffuse faster due to the strong dependence of the available volume factor Vm,n on to + n, which means that most of 
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Jt(^) -^ ^ ^ ^ 

Jb{x) —^ ^ ^ . . /3 = l/(5min) 

Jw{x) —^ ^ ^— ^ <-•-•■ 


- -t ^ ^ -t -• 

-i -1 -I • - - ■ ' /3 = 1/(10 min) 

H «-<-<-<- {- I- 


. /3 = l/(20min) 


0 0.25 0.5 0.75 1 

x/l 

FIG. S3: Fluxes of translating (T), transiently bound (B) and free (F) ribosomes along the cell’s long axis: The arrow length is 
proportional to local ribosome flux, and the arrows in the legend correspond to a flux of 20/s. The fluxes are shown for different 
values of mRNA degradation rate j5, where j3 and the total mRNA production rate Qtot are varied together, keeping the total 
mRNA number constant. 


the diffusive flux of mRNAs is due to polysomes with a number of ribosomes substantially smaller than average. Thus 
our use of the average ribosome loading (IVt + iVB)/fVmRNA in Eq. (S48) substantially overestimates the ribosome flux 
carried by mRNAs. However, the estimate J^+b in Eq. (S48) does correctly capture the full numerical result that 
the ribosome flux is proportional to the total mRNA production rate atotj see Fig. S3. This proportionality leads us 
to conclude that ribosome circulation is driven by the flux of new mRNAs from the nucleoid to the poles, and not 
by binding of mature mRNAs and F ribosomes in the nucleoid with subsequent expulsion from the nucleoid due to 
excluded-volume effects as discussed in the main text. 


S9 Model including nucleoid-bound mRNAs 


In this section, we extend our model to include the fact that mRNAs remain bound to the nucleoid while being 
transcribed. We introduce an additional set of nucleoid-bound mRNA species with m B ribosomes and n T ribosomes, 
whose densities we denote by p’^^{x,t). The reaction-diffusion equations for the nucleoid-bound mRNA densities are 


- konCF{x,t)p*^ „^{x,t) - k^smp*^ ,^{x,t) - kl^c^{x,t)p*^ ,^{x,t) - kl^n p*^^^{x,t) 

+ fconCF(a:, t) + k^^ {m + 1 ) p*m+i,nix, t) + kl^c^{x, t)p*^^„_^{x, t) + k'^ff {n + 1 ) t) 

+5m,o5n^Qa{x) - Pp^ „{x,t) -'yp*^^^{x,t). (S49) 

Since the mRNA species p*^ „ are bound to the nucleoid, there is no diffusion term. The terms in the first and second 
lines on the RHS correspond to ribosome binding and unbinding as in Eq. (2). The first term in the last line represents 
mRNA transcription, the second term co-transcriptional degradation [13], while the last term represents unbinding of 
an mRNA from the nucleoid at rate 7 . The dynamics of free mRNA species, i.e. mRNAs that are not bound to the 
nucleoid, is still described by Eq. (2), where the niRNA-transcription term in the last line is replaced by the source 
term TpJ/j which represents a nucleoid-bound mRNA becoming a free mRNA: 


dptn,n{x,t) 

dt 


dpm,n{x,t) 

dt 


D 


d'^Pm,n{x,t) 

dx‘^ 




d'^Vm+n{x) 


dx^ 


-k^^C^{x,t)pm,n{x,t) - k^gmpm,nix,t) - kl^C^{x,t)Pm,n{x,t) - k'^g n Pm,n{x,t) 

T (^7 f )Pm—l,n f^off T ^on^F (^5 1 (^; 0 T ^off T 0 

+1 - P Pm.ni^^t)- (S50) 


We take the unbinding rate of nucleoid-bound mRNAs to be 7 = 1/min, i.e. the inverse of the average time for mRNA 
transcription, compare section II. Finally, the dynamics of F ribosomes is described by the analog of Eq. (1): 
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dcF{x,t) 

Wt 


= D 


d'^Vj.{x) 

-Vp{x) — Cp{x,t)- 


dx^ ^^2 

^max 1 ^max 

E E [Pm,«(a;, t) + p*ra,nix, t)] + k^s E E 


^^max "^max 


m—0 n—0 

^max ‘^max 1 


m—1 n =0 

^^max "^max 


E E [Pm,«(a;, t) + p^,„(a;, t)] + klg E E n[p^^ri{x,t) + p*^^^{x,t)] 

m—0 n—0 m—0 n—1 

^^max ^max 

+/3 EE (TO + n)[p,„,„(x,i) + Pm,„(a;,t)]. (S51) 


m—0 n—0 


We now consider the steady state of the equations above for Pm,n, Pm n and we assume that B ribosomes 

are in rapid equilibrium with F ribosomes. Proceeding along the lines of section S 6 , in the rapid-equilibrium limit we 
have 


Pm,n\^) 

Pm,n(.x) 


1 

ml 

1 

ml 


( feonCF(x) 

V fcoff 


( konCFjx) 
\ ^oS 


Pn{x), 

Pnix), 


(552) 

(553) 


where Pn(x) = po,n(a^), and Pn{x) = Pq„(x). In the rapid-equilibrium limit, the densities Pn(x) satisfy the following 
reaction-diffusion equation 


-I- 


f kon^F^x) , . , .. 

J ^Fla:)] 

dpn{x) 


^Pn{xW{x) 


dcp(x) , , , .dVf{x) 

Vp{x) — Cf(X) 


dx 


/ \ / \dv^{^x^ , . , , 

Vn{x) - Pn[x)— - n Pn{x)v^[x)- 


dx 

1 dvp{x) 


(S54) 


-I- exp 


konCFjx) 

r.B 


dx dx dx 

{-konCF{x)pn{x) - k'^g uPn{x) + k’^^c^{x)Pr,-i{x) + kjg {u + 1 ) p„+i(a;) - Ppn{x) + 'yPnix)] = 0 , 


while Pnix) satisfies the following equation involving only reaction terms 
kon^Fix) 


- j [-KnCFix)p*nix) - k'^gn p*nix) + k'^nCFix)Pn_iix) + k^g in + 1) Png^ix) - iP + j)Pnix)] (S55) 

+aix) = 0 . 


Finally, the rapid-equilibrium reaction-diffusion equation for Cf(x) reads 


+^;£”E j exp ) [up(a:)]^ 


n—0 


^off 


i-f 


fc^„Cp(x) 

'^off 




Up (a;) 


d^cJx) , . , .d?vJx) 

, 2 W(a^) - Cp(x) 

dx‘^ dx^ 


dCfix) 

dx 


Vfix) — Cp(x) 


dvpix) 


k^nCFix) (dpnix) dvnix)\ ^?nCF(a;) , . , . 

^ 7 B H” , 3 ( J ^r(^) Pn\^) 7 I Vf\^) , g Pn\^)'^K\^)'^ 

Kg Kg \ dx dx J k^g 


exp 


k^nCFixY 


dx 

^max 1 


dx 

dvpix) 

dx 


uB 

'^oG 


n—0 


- konCFix) + P*Kx)] + fcoff E + P*nix)\ 

n—1 

^max I 

+/3 E '<^\pKx) + P*nix)] \ 


(S56) 


= 0 . 


n—0 


For both free and nucleoid-bound mRNAs, we enforce no-flux boundary conditions at the cell pole and at the cell 
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center 




dx 

dp*n{x) 

dx 


WR(a;)[t>p(x)]” -p*(x) 


dx 

dx 


= 0 , 


x—0,£ 


= 0 . 


(557) 

(558) 


x=0,i 


For ribosome concentrations, we impose the no-flux boundary condition at the cell pole, Eq. (S39), and a constraint 
on the total number of ribosomes 

o£ C CXD Tlmax ^ 

Ntot = 2 dx <c^{x) + ^ ^{m + n)[p™,„(a;) -f p^,„(a;)] > 

do { m=0 n=0 ) 

= 2 dx |cp(x) + exp + n) [p„(x) + ^^(x)]| . 


Z-B I ^ \ 

/ n=0 ^ '‘‘off 


(S59) 


Finally, the ribosome efflciency reads 


E = 


OO Tlmax p£ 

2kis E E / dx[pm,n{x) + P*m,ni.x)] 

^—n ^ — 1 0 


m=0 n—1 


— 2 n dx exp 


konCv{x) 


[Pn{x) + pI{x)]. 


(S60) 


We solved numerically Eqs. (S54), (S55), and (S56) for pn, p* and Cp, by fixing the maximal number of allowed T 
ribosomes per mRNA at Umax = 24. The resulting mRNA profiles and ribosome concentrations are shown in Figs. S4 
and S5, respectively. The results confirm the qualitative behavior of the model with only free mRNAs, see Figs. 3 and 
4. The plots of ptot (a^) and p^^^ ix) = show that the large majority, i.e. ~ 83%, of 

mRNAs are free, while the remaining ~ 17 % of mRNAs are bound to the nucleoid, where the latter fraction constitutes 
also an estimate for the percentage of mRNAs that are degraded co-transcriptionally [13]. Free mRNAs are strongly 
segregated from the nucleoid, while nucleoid-bound mRNAs, are by construction localized in the nucleoid region. For 
free mRNAs, the larger m, n, the stronger the segregation away from the nucleoid, and the overall free polysome 
distribution is still peaked around ~ 10 T ribosomes and ~2 B ribosomes per mRNA as in Fig. 3. For nucleoid-bound 
mRNAs, the vast majority of mRNAs have low m, n, with the distribution peaked at m ~ 1 and n ~ 3 — 4, reflecting 
the somewhat lower concentration of F ribosomes in the nucleoid (Fig. S5) and, more importantly, the limited time 
for T ribosomes to bind before the mRNAs become free. The localization pattern of these typical bound mRNAs 
matches the profile of the nucleoid, i.e. the profile of mRNA production. By contrast, the spatial distributions of the 
few bound mRNA species with large m, n are strongly peaked at the edge of the nucleoid, reflecting the higher density 
of free ribosomes outside versus inside the nucleoid. 

Figure S5 shows the steady-state concentrations, Ct(x) and Cb(x), of T and B ribosomes loaded on free mRNAs, 
given by Eqs. (S42) and (S43), the concentration of F ribosomes Cf(x), as well as the fluxes Jt(x), Jb(x) given by Eqs. 
(S45), (S46), and Jf(x). The figure also shows the concentrations of T and B ribosomes loaded on nucleoid-bound 
mRNAs, c* (x) = E“=o Y.n=TP*m,u(^) and c* (x) = Em=i E”™ respectively. While T and B ribosomes 

loaded on free mRNAs are strongly excluded from the nucleoid, F ribosomes can penetrate the nucleoid region. In 
addition, the concentration profiles above provide an estimate of the fractions of ribosomes in the nucleoid region that 
translate mRNAs co-transcriptionally and post-transcriptionally. The number of ribosomes that translate co- and post- 

transcriptionally is Wo = 2 Em=o En”!" ^ dxp^_„(x) ~ 3.6 x 10^ and iVpost = 2 Em=o E””!" dxpm,n{x) ~ 

3.9 X 10^, respectively, while the total number of ribosomes in the nucleoid region Wn ~ 1-1 x lO"* can be obtained 
by replacing £ with £/2 in the expression (S59) for the total number of ribosomes in the cell. Thus, Wo/Wn ~ 34% 
of ribosomes carry out co-transcriptional translation in the DNA-rich region, while a fraction of Wost/-^in ~ 37% 
translates post-transcriptionally. As shown in Fig. S4A, the larger the loading number of a free mRNA, the stronger 
its nucleoid segregation: it follows that such post-transcriptional translation in the nucleoid occurs mostly on free 
mRNAs loaded with a relatively small number of ribosomes, which could be either newly transcribed mRNAs that 
did not have time to diffuse out of the nucleoid, or old mRNAs with small loading number which penetrated into the 
nucleoid. Figure S5 also shows that there is a net poleward flux of T and B ribosomes loaded on free mRNAs and a 
net flux of F ribosomes toward the cell center, in agreement with the model with only free mRNAs. Finally, T and B 
ribosomes loaded on nucleoid-bound mRNAs are localized approximately uniformly throughout the nucleoid region. 
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FIG. S4: Steady-state mRNA and polysome distributions for the extended model including both free and nucleoid-bound 
mRNAs, where we hxed the maximal number of allowed translating (T) ribosomes per mRNA at n < rimax = 24, and we show 
mRNA species with m < 8 transiently bound (B) ribosomes. (A) mRNA and polysome distribution for free mRNAs. Total 
mRNA density (red) and density pm,n{x) of mRNAs with m B ribosomes and n T ribosomes (gray). The density po,o(x) 

of ribosome-free mRNAs (green) and the density P 8 , 24 (x) of mRNAs with the largest number of B, T ribosomes considered 
(blue) are also shown. The profiles pm,n{x), po,o(x}, and P 8 , 24 (x) are normalized to unit area. Inset: distribution of mRNA 
species, shown as a heat map of the number Nm,n = /g dxpm,n{x) of mRNAs with m B ribosomes and n T ribosomes in the 
right half of the cell. The maximal number of allowed T ribosomes per mRNA, rimax = 24, is also marked. (B) Same as panel 
A for nucleoid-bound mRNAs, with „ dxp^ „(x) . 


A 



Jt(x) —^ 

Jb{x) —^ 
Je(x) 






x/£ 


FIG. S5: Steady-state ribosome concentrations and ribosome fluxes for the extended model including both free and nucleoid- 
bound mRNAs. (A) Goncentrations and fluxes for translating (T) and transiently bound (B) ribosomes loaded on free mRNAs, 
and for free (F) ribosomes. On the top, we show the concentrations Ct{x), Cb(x), and CF(a:) of T, B, and F ribosomes, for 
the right half of the cell. Bottom: Fluxes of T, B, and F ribosomes along the cell’s long axis depicted in the top panel, 
where the arrow length is proportional to local ribosome flux, and the arrows in the legends correspond to a flux of 30/s. (B) 
Goncentrations c^(x), Cb(x) of T and B ribosomes loaded on nucleoid-bound mRNAs. 


We next make use of the model with nucleoid-bound mRNAs to study the contribution of co-transcriptional trans¬ 
lation to ribosome efficiency. To do so, we remove co-transcriptional translation by fiat in our model by eliminating 
binding of T ribosomes to nucleoid-bound mRNAs, but still allowing rapid-equilibrium binding and unbinding of B 
ribosomes to these mRNAs. First, without co-transcriptional translation the first four terms on the left-hand side 
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(LHS) of Eq. (S55) vanish: as a result, one has 

P*n{x) = 5„,oa(^)/(/3 + 7). (S61) 

Substituting this in the reaction-diffusion equation (S54) for Pn{x), we recover the intuitive result that the reaction- 
diffusion equations and boundary conditions for the two mRNA species with no co-transcriptional translation are 
equivalent to those for free mRNAs only, Eq. (S27), with a reduced rate of mRNA production, a{x) —)■ a{x)^/ (/3 -I- 7 ), 
accounting for the decay of some mRNAs before transcription is complete. Second, in the absence of co-transcriptional 
translation the terms proportional to p* in the reaction-diffusion equation (S56) vanish, and the reaction-diffusion 
equation for Cp also reduces to Eq. (S37) for the model with free mRNAs only. Finally, let us consider the boundary 
conditions of Eq. (S56) in the absence of co-transcriptional translation: While the no-ffux condition at the cell pole, 
Eq. (S39), still holds, the constraint (S59) on the total number of ribosomes should be modified to take account of 
the fact that there are no T ribosomes on nucleoid-bound mRNAs: 
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(S62) 


where in the second line we used Eqs. (S52) and (S53), and in the third line we substituted Eq. (S61). The constraint 
(S62) for the total ribosome number and its analog (S40) for the model with no nucleoid-bound mRNAs differ by the 
term proportional to a{x)'. intuitively, this term represents the total number of B ribosomes carried by mRNAs that 
are being transcribed in the nucleoid. 

We will now estimate the ribosome efficiency, thus quantifying the change in efficiency due to the loss of co- 
transcriptional translation. Without co-transcriptional translation, ribosomes are not allowed to translate nucleoid- 
bound mRNAs, thus the expression for the total rate of protein production S is given by Eq. (S47), and the ribosome 
efficiency is given by e = S/Wot = 1-88 x 10“^/s, which is ~3% smaller than the efficiency e = 1.94 x 10“^/s in the 
case with co-transcriptional translation. This loss of efficiency is mostly due to the fact that nucleoid-bound mRNAs 
carry B ribosomes: indeed, compared to the case with co-transcriptional translation, this additional fraction of B 
ribosomes implies a reduced pool of F ribosomes, thus decreasing the F —)■ T rate, and consequently the ribosome 
efficiency. 


SIO Model including SOS and 50S ribosomal subunits 

In this section, we extend our model to include the fact that each bacterial ribosome is composed of a 30S and a 50S 
subunit. In section SIO A we introduce the reaction-diffusion equations for the two-subunit model, whose parameters 
will be discussed in section SIOB. In section SIOC we discuss the rapid-equilibrium approximation for B 30S and 50S 
subunits, and in section SIOD we present the numerical solution of the resulting reaction-diffusion equations. 


A. Reaction-diffusion equations 


The reaction-diffusion equations for the mRNA densities are 
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In both Eqs. (S63) and (S64), the terms in the first line represent mRNA diffusion, where vi+m+n, r’;_|_^_i_„ are 
the available-volume profiles for mRNAs without or with a 30S subunit occupying the translation-initiation site, 
respectively. The second, third, fourth, and fifth line represent binding and unbinding of 30S and 50S subunits at 
rates k^^ and fc®®®, /c^fj® respectively, compare Eq. (2). The terms in the sixth line represent unbinding of 

70S pairs [14] at rate A:^^®. The first two terms in the seventh line represent a 30S subunit binding and unbinding 
the mRNA translation-initiation site at rates k'A^^, k'J^^ respectively, while the third term represents a 50S subunit 
binding to the 30S subunit on the initiation site at rate k'FP^ and forming a T 70S pair. Finally, the terms in the last 
line of Eq. (S63) represent mRNA transcription and degradation, and the term in the last line of Eq. (S64) represents 
mRNA degradation. The reaction-diffusion equation for the concentration of free 30S subunits reads 
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l—O m—0 n—0 


In Eq. (S65), only mRNA species with allowed values of 0 < / < Imaxi 0 < m < rUmax) and 0 < n < rimax are 
considered. The terms in the second line represent 30S subunits transiently binding and unbinding mRNAs, compare 
Eq. (1). The third line represents unbinding of a 30S subunit coming from a 70S pair, while the two terms in 
the fourth line represent a 30S subunit binding and unbinding the initiation site, respectively. Finally, the last line 
represent 30S subunits being freed from mRNA molecules as these are degraded. The reaction-diffusion equation for 
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the concentration of free 50S subunits reads 
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The terms in the second line represent SOS subunits transiently binding and unbinding mRNAs. In the third line, the 
first term represents unbinding of a SOS subunit coming from a 70S pair, and the second term represents a SOS subunit 
binding the SOS subunit on the initiation site. Finally, the term in the fourth line represents SOS subunits being freed 
from degraded mRNAs. 

We consider Eqs. (S63)-(S66) at steady state, and in Eqs. (S63), (S64) we set a no-flux boundary condition at the 
cell pole and at the cell center: 
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In Eqs. (S6S), (S66), we set no-flux boundary conditions at the cell pole 
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and we enforce a constraint on the total number of SOS and SOS ribosomal subunits 
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B. Model parameters 


The model parameters in Eqs. (S63)-(S72) are fixed as follows: The diffusion coefficients D, Dp, the degradation 
rate /3, and the ribosomal available volume Vp{x) are taken to be the same as in the model with one ribosomal unit, 
see section II. Assuming for simplicity a radius of 10 nm for both the SOS and SOS subunits, the available volume 
vi+rn+nix) is given by Eq. (S12). It follows that the available volume for mRNAs with an additional SOS subunit at 
the initiation site is vAm+nA) ~ vi+m+n+i{x). We fixed Aaostot and A^sostot from the total number of ribosomes 
and from the 30S:S0S stoichiometric ratio, which can be estimated as follows: First, every SOS or SOS subunit includes 
one 16S or one 23S ribosomal RNA (rRNA), respectively. As a result, the 30S:S0S ratio can be estimated from the 
stoichiometry of 16S- and 23S-rRNA transcript abundances, which is found to be close to 1:1 [IS]. Second, the 1:1 
ratio above is generally supported by studies of the abundance of ribosomal proteins in the two subunits. Third, the 
experimental observation that ~ 80% of SOS ribosomal subunits are engaged in translation [S], with similar results 
for SOS subunits [16], implies that the great majority of ribosomes must come as pairs of SOS and SOS subunits, thus 
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confirming a 30S:50S ratio close to 1:1. We thus assume a 1:1 ratio of 30S:50S subunits. Consequently, for a total 
ribosome number of Wot = 6 x 10^ [5], we obtain Wostot = Wostot = Wot- The transient-binding rates 
^on^i ^off^ t)e assumed to be significantly faster than any other relevant time scale, and they will be treated in 
the rapid-equilibrium limit as in section S6. In section SIOC, we will show that in this limit the reaction-diffusion 
equations depend only on the ratio between the binding and the unbinding rates, which can be computed along the 
lines of section II, providing the estimate ~ ^oS^/^off^ — 5.4 x 10“^ /rm. The binding constant of 30S 

subunits on the initiation site is estimated as follows: An in vitro kinetic analysis of the assembly of the 30S initiation 
complex quantified the binding rate of 30S subunits to a single mRNA as fcan c^o, i.e. the rate is proportional to the 
three-dimensional 30S concentration Can, and to an associated rate constant kso- The value of k^^ was estimated from 
the slope of the 30S-subunit binding rate as a function of the three-dimensional 30S concentration, and was found to 
be k^r) ~ 0.1 /xm^/s [17]. The 30S binding rate in our one-dimensional model is then obtained by dividing fcao by the 
area of the two-dimensional cell slice perpendicular to the cell axis: k'J^^ = k^Y>/= 0.13 /xm/s, where R = 0.5 /xm 
is the radius of a circular cell slice. The unbinding rate of 30S subunits from the initiation site A^l^® = 2/s is taken 
from the same in vitro kinetic analysis [17]. The rate of unbinding of 70S pairs can be estimated along the same lines 
as /cjfj in the model for one ribosomal unit in section II, yielding /c/jj® = kjg. Finally, the binding rate of 50S subunits 
to 30S subunits on the initiation site can be obtained from the other parameters by imposing the conservation of the 
total number of 30S and 50S subunits: 

klfN,,s - fc/f ®Wos pP + k'J°^£7 = 0, (S73) 

klfN,os-k',l°^N,os.p' = 0, (S74) 

where Wosfj Wosfi and Wos are the total number of F 30S and 50S subunits and 70S pairs, and p, p' are the average 
total axial densities of mRNAs without or with a 30S subunit on the initiation site, respectively. Equation (S73) 
imposes the balance between the overall binding and unbinding rates of 30S subunits in the cell: The first term in the 
LHS represents the total rate at which 30S subunits in a 70S pair unbind from mRNAs, while the second and the third 
term are the total rates at which 30S subunits bind and unbind from the translation-initiation site, respectively. Note 
that the rates of binding and unbinding of transiently bound 30S subunits do not appear in Eq. (S73) because the sum 
of these rates is zero. Denoting by Wos' the number of 30S subunits bound to the translation-initiation site, in Eq. 
(S73) we omitted the rate /3(Wos + Wos') at which 70S subunits and 30S subunits bound to the translation-initiation 
site are freed by mRNAs that are degraded, because this rate is significantly smaller than kl^N^gs- Equation (S74) 
represents the balance between the total binding and unbinding rates of 50S subunits. Similarly to Eq. (S73), we 
omitted the term /3Wos 'C A:/[j®Wos associated with mRNA decay. The rate fc/^°® at which 50S subunits bind to 
30S subunits on the translation-initiation site is then determined as follows: We set Wos f = Wos f = 2 % Wot and 
Wos = 80% Wot [5, 16], we use the fact that the mRNA densities p and p' are related by the constraint on the total 
mRNA number p + p' = A^mRNA/(2^), and we solve Eqs. (S73) and (S74) for p, p', and fc/^®. 


C. Rapid equilibrium for transiently bound SOS and SOS subunits 


In this section we will derive the reaction-diffusion equations assuming that B 30S and 50S subunits are in rapid 
equilibrium with their F counterparts. We present only the main steps of the derivation, which follows along the lines 
of the rapid-equilibrium approximation for the model with one ribosomal subunit discussed in section S6. 

In the rapid-equilibrium limit, we set 
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where A ^ 1 and the boldface binding-unbinding rates are of order unity. 

To leading order in A, we obtain 

[fcoff®(l -k l)pi + i^m,nix) - kl°^C30SFix)pi,m,n{x)] - I ppm,nix) “ ^on® C 30 S f (a:)pi-l,m.n (a^)] 

T[^off (jXl l)p/^yj3+l,n(a:) ^on ^503 F(a:)p/^m,n(a:)] [^ofF xn pp'm,n(yX^ ^on C 50 S F 1-''^ (^)] 


(575) 

(576) 


0. (S77) 


The four terms in the LHS of Eq. (S77) represent the mRNA current between species (?,to) and the four species 
{I + 1,to), {I — l,m), {l,m + 1) and {l,m — 1), thus implying that the sum of the four currents is zero. Note that this 
condition does not necessarily imply that each current vanishes. Indeed, there could be solutions pprn,n{x) involving net 
nonequilibrium currents between pairs of mRNA species, that would still satisfy the current-balance condition (S77). 
However, if we assume that the 30S and 50S transient binding and unbinding is a purely spontaneous, equilibrium 
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process that is not coupled to any external drive or energy source, all the individual currents between mRNA species 
must vanish 


+ + - kl°^C30SF{x)pi,m,n{^) = 0, ^<^max-l, (S78) 

+- A:on®C5osF(a;)p;,m,n(a;) = 0, m < m^ax - 1- (S79) 

Equations (S78) and (S79) can now be solved along the lines of Eq. (S19) for one ribosomal unit, and their solution 
is analogous to Eq. (S20), providing the following leading-order expression for pi^rn,n{x): 
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Proceeding along the same lines, we obtain the leading-order expression for pj ^ n(^)- 
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Next, we derive the reaction-diffusion equations for 
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Proceeding along the lines of Eqs. (S21)-(S27), we introduce the local fluxes of mRNAs 
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We sum Eqs. (S63) and (S64) with respect to ^ = 0,..., /max and m = 0,..., rumax, and we obtain the following 
reaction-diffusion equations to leading order in A 
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Proceeding along the lines of Eqs. (S26) and (S27), we rewrite Eq. (S88) as a reaction-diffusion equations for p„(x): 
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Similarly, Eq. (S89) implies the following reaction-diffusion equations for Pn{x): 
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dx 


' exp 


f kon^C30SF{x) k^^^C30SF{x)\ r 7,703^ J ( 

^ ^303 i-sns ) L ^off ^Pn\X) 


7,503 
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WF(a;) dx 

+kls^{n + l)p'„+i(x) -f ®C3osF(a;)p„(a;) - k'Ji°^p'^{x) - k'J,°^C3osF{x)p'r,{x) - I3p'^{x)] = 0. 

We now derive the rapid-equilibrium reaction-diffusion equation for the ribosome concentration C 3 osF(a^): We intro¬ 
duce the subleading corrections to pprn,n and p'^ ^ setting 


Pi 


,m,nix') — Ppm,nix^ -|- ^pi^jyi^nix^ 


Pi,m,nix) = -f -Ap; „(a;). 


(592) 

(593) 


We rewrite the steady-state reaction-diffusion equations (S63) for pi,m,nix) as follows: 

[feoff® {I + 1 ) Api + I^rn,nix) - feo°®C 30 S f( a^) App,„_„ (x)] - [fe^ff® I Api^rn,nix) - feo°®C 3 os F(a;) Ap/_i,„,„(a;)] 
Pl,m-\-l,nix) ^on ^50SFix)pi,ni,nix)] [^off ^ Pl,ni,nix) ^on l^50S pix)Pi,ni—l,nix)] 
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kon ^iospix)Pl,m,nix) ~\~ feofT Pl,m,nix) “f ^on ^50S F X (a:) -f dl^Q dm ,0 dn,Q Ol(x) [d Pl,m,nix) 


0, (S94) 


where in the first two lines we collected all the terms involving fast processes related to B 30S and 50S subunits, and 
in the last two lines we collected all the remaining terms describing slow processes, such as diffusion and unbinding of 
70S pairs. Also, in the first line we used Eq. (S92), and we observed that the terms to leading order in A cancel out 
because of Eq. (S80). To derive the reaction-diffusion for C 3 osF(a^), we sum Eq. (S94) with respect to m = 0,..., rumax, 
and we obtain to leading order in A: 


[klfl Apf°„®(x) - fe^f c3osF(a:)Ap?°®,„(x)] - [kl^i{l + l)Apf°l„(x) - fe^f C3osF(a:)Apf°®(x)] = 
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+ 5l^Q5n,OOl{x), 


(S95) 


where 


"^max 

^Ppnix) = E ^Pk^Ax), (S96) 

m—0 

and the superscript ‘30S’ refers to the fact that Apf^^ is the total subleading density of mRNAs with I 30S subunits, 
i.e. it results from a sum over all contributions of 50S subunits in the RHS of Eq. (S96). Intuitively, by summing Eq. 
(S94) with respect to m, we averaged out the contribution of SOS subunits, thus obtaining the overall current between 
mRNA species with I 30S subunits and species with I — 1 and I + 1 30S subunits, compare the two terms in the LHS 
of Eq. (S95), respectively. Given that the contribution of SOS subunits has been averaged out, Eq. (S9S) now has 
the same form as Eq. (S30) for a single ribosomal unit, and it can thus be solved iteratively along the lines of Eqs. 
(S31)-(S32): the result is 

^max '^max 

klflApf^^ix) - fe3°®C3osF(a;)Apf°®i „(a;) = ^ ^ 

p—l m—0 

Importantly, to derive the reaction-diffusion equation for C 3 osF(a:), it is not necessary to determine all the currents 
between mRNA species (?,m) and species (^ + l,m), (^—l,m), {l,m + l) and — that appear in Eq. (S77): the 
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(S97) 
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average current in Eq. (S97) is all we need. Indeed, let us consider the second line of Eq. (S65): the terms involving 
Pi,m,n can be rewritten in terms of the current in Eq. (S97) as follows: 
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(S98) 


where in the second line we used Eqs. (S80) and (S92) and we observed that the terms of order A cancel out, in the 
third line we used Eq. (S96), in the fourth line we used Eq. (S97), and in the fifth line we used the identity (S34). 
Similarly, the terms in the second line of Eq. (S65) involving pj ^ „ read: 
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(S99) 


We now compute explicitly the RHS of Eqs. (S98) and (S99) for large ^maxi along the lines of Eqs. (S35) and 

(S36), we substitute the result in Eq. (S65), and we obtain the final reaction-diffusion equation for the concentration 
of F 30S subunits 
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In Eq. (SlOO), the terms in the second, third, and fourth lines represent the diffusive flux of B 30S subunits carried 
by diffusing mRNAs with no 30S subunits on the initiation site, and the terms in the fifth, sixth and seventh lines 
represent the diffusive flux of B 30S subunits carried by mRNAs with a 30S unit at the initiation site. 
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The reaction-diffusion equation for the concentration of F 50S subunits can be obtained along the same lines as Eq. 
(SlOO), and it reads 
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The boundary conditions for Eqs. (S90), (S91), (SlOO) and (SlOl) are obtained from Eqs. (S67)-(S72) and from 
the leading-order expressions (S80) and (S81) for the mRNA densities. The boundary conditions for Eqs. (S90) and 
(S91) read 
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The no-ffux boundary conditions for Eqs. (SlOO) and (SlOl) are given by Eqs. (S69) and (S70), while the constraints 
(S71), (S72) on the total number of SOS and 50S ribosomal subunits read 
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Overall, Eqs. (S69), (S70), (S90), (S91), (SlOO), (SlOl) and (S102)-(S105) characterize the steady-state behavior of 
the reaction-diffusion model with two ribosomal subunits in the rapid-equilibrium limit. Compared to the model with 
a single ribosomal subunit discussed in section II, this two-submit model provides a more realistic characterization of 
the transcriptional-translational process, describing explicitly the assembly and disassembly of 70S pairs on mRNA 
molecules. More realistic models describing the internal degrees of freedom of ribosomes and mRNAs could be obtained 
with full molecular-dynamics simulations, rather than with reaction-diffusion equations for local concentrations. Even 
if they are likely to be limited to a single ribosome and mRNA molecule, such molecular-dynamics simulations may be 
useful for describing several interesting features of the transcriptional-translational machinery on a microscopic level, 
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such as recycling of ribosomes which translate mRNA molecules forming loops, where initiation and termination sites 
are in close proximity [18]. 

Proceeding along the lines of section S 6 , all physical quantities of interest can be expressed in terms of p„, C 30 SF) 
CsosF by using Eqs. (S80)-(S83). Some of these quantities are the total mRNA densities 
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and the concentration of 30S subunits bound to the translation-initiation site 0333 ( 0 ;) = p^^^(x). Other quantities are 
the fluxes of F 30S and SOS subunits, that we denote by Jsospix) and o/sosf)^;), and that are obtained from Eq. (S3) 
by replacing Cp with C 30 SF and Csosfi respectively, the flux of 70S subunits 
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and the flux of SOS subunits bound to the translation-initiation site 
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D. Results 


We numerically solved Eqs. (S90), (S91), (SlOO), and (SlOl) by fixing the maximal number of allowed translating 
70S pairs per mRNA at nmax = 24, and the results are shown in Figs. S, S18, and S7. In addition to the features 
of these results discussed in section II, here we observe that the total mRNA distributions in Fig. S18 have a dip at 
xj^ ~ 0.4. This minimum is due to the fact that the total mRNA profile is the superposition of two peaks: One peak in 
the nucleoid region due to mRNAs that are being transcribed and are nearly free of ribosomal subunits, and one peak 
at the cell poles due to mRNAs loaded with multiple subunits, which are strongly excluded from the nucleoid. In Fig. 
S18, there is an average number of ^ 14 subunits, i.e. 30S, 50S and 70S pairs, per mRNA, implying excluded-volume 
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FIG. S6: Steady-state mRNA distributions for the model including SOS and SOS ribosomal subunits, where we fixed the maximal 
number of allowed translating (T) 30S-50S (70S) pairs per mRNA at n < rimax = 24, and we show mRNA species with I < 4 and 
m < 4 transiently bound (B) SOS and SOS subunits, respectively. (A) mRNA profiles for mRNAs with no SOS subunit bound 
to the translation-initiation site. Total mRNA density ptot{x) (red) and density pi,m,n{x) of mRNAs with I B SOS subunits, 
m B SOS subunits, and n 70S pairs (gray), where we show only the curves with I, m and n even for greater clarity. The 
density po,o,o(x) of ribosome-free mRNAs (green) and the density p 4 , 4 , 24 (a:) of mRNAs with the largest number of B, T subunits 
considered (blue) are also shown. The profiles pi^m,n{x), po,o,o(x), and P 4 , 4 , 24 (x) are normalized to unit area. (B) Same as panel 
A for mRNAs with a SOS subunit bound to the initiation site. 
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FIG. S7: Polysome distributions for the model including SOS and SOS ribosomal subunits, where we fixed the maximal number 
of allowed translating (T) S0S-50S (70S) pairs per mRNA at n < rimax = 24, and we show mRNA species with I < 4 and m < 4 
transiently bound (B) SOS and SOS subunits, respectively. (A) distribution of mRNA species with no SOS subunit bound to the 
translation-initiation site, shown as a heat map of the number Ni^rn,n = Jq dxpi^m,n{x) of mRNAs with I B SOS subunits, m 
B SOS subunits, and n T 70S pairs in the right half of the cell. The diameter of each dot is proportional to Ni^rn,n, where the 
largest dot in the plot corresponds to Ni^rn,n ~ 1- (B) Same as panel A for the number = fo dxp'i „^ „(x) of mRNAs with 

a SOS subunit bound to the initiation site, where the largest dot in the plot corresponds to „ ~ 20. 

effects strong enough that the second peak shrinks toward the cell poles, thus leaving a dip between the cell center 
and the poles—see also section S14. 

We estimate the total translation rate in our model with two subunits as 

5 : 5 : 71 / dx[ppmAx) + P'l,m,n{x)] (S115) 

/=0 m=0 n=l "^0 

and the translation efficiency as e = S/A^sostot = 1-94 x 10“^/s, which is practically identical to the efficiency 
e = 1.94 X 10“^/s for the model with a single ribosomal unit discussed in section 11. 

We conclude this section by discussing how our results pertain to in vivo conditions. In this regard, we recall 
that all the model parameters discussed in section SlOB have been estimated from in vivo observations, except the 
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FIG. S8: Translation-rate increases d'E/dN^ostot and dE/dN^ostot due to the addition of a SOS and a 508 subunit, respectively, 
as functions of the unbinding rate of SOS subunits from the translation-initiation site. The in vitro value of the unbinding 

rate = 2/s is marked on the top horizontal axis. 


binding-unbinding rates of SOS subunits to the translation-initiation site, and which have been estimated 

in vitro [17]. First, we note that the higher concentration of solvated molecules present in vivo may significantly 
affect the binding-unbinding rates: a relevant example is that of DNA-bound proteins, whose dimeric nature allows 
for partial dissociation and thus enhanced exchange in the presence of competing proteins in solution [19-21]. It 
is natural to hypothesize that the same mechanism might apply to SOS subunits bound to the translation-initiation 
site, resulting in a significant increase in the unbinding rate in vivo compared to the in vitro value k'Jg^ = 2/s 

discussed in section SIOB. Second, using the in vitro values of k'J^^ and and the solution of Eq. (S73), (S74), 
we obtain that the fraction of mRNAs with a SOS subunit bound to the translation-initiation site is p'/{p + p') ~ 96 %, 
i.e. the vast majority of mRNAs. It follows that the overall translation rate E is mostly limited by the number of 
SOS subunits rather than by the number of SOS subunits. Indeed, let us consider the quantities dE/dN^ostot and 
dE/dN^ostot, which represent the translation-rate increase due to the addition of a SOS or a SOS subunit, respectively. 
These two translation-rate increases are shown in Fig. S8 as functions of the unbinding rate For values of 

the unbinding rate close to the in vitro value, we have dE/dN^ostot ^ dE/dN^ostot- in this case, there is a “waste” 
of SOS subunits, i.e. the number of SOS subunits could be significantly reduced without significantly affecting the 
translation rate. On the other hand, for larger values of > 10^/s, the two translation-rate increases are similar, 
i.e. i9 E/(97V3 os tot ~ 5E/i9A^5os tot : in this case, the SOS and SOS subunits are co-rate limiting for translation, reflecting 
an efficient allocation of cellular resources which is generally expected in vivo [22, 2S[. Overall, these observations 
provide an example of how our model can be used to relate microscopic parameters such as the binding and unbinding 
rates of ribosomal subunits to macroscopic observables such as subunit stoichiometries and ultimately to functional 
performance, namely the protein-translation rate. 


Sll mRNA degradation by RNAse enzymes 


In this section, we incorporate in the two-subunit model above a more realistic description of mRNA degradation 
based on the mechanism of mRNA decay suggested by in vivo studies in E. coli: mRNA degradation is triggered by 
RNase enzymes, which target the 5' end of mRNA and compete with ribosomes to bind the transcript and degrade it 
[1S[. In addition, RNase enzymes have been suggested to localize to the cell membrane [24]: In our one-dimensional 
reaction-diffusion model, such a distribution of RNase on the three-dimensional membrane will result in an effective, 
one-dimensional distribution of RNase on the long cell axis. To obtain this distribution, we model the cell as a cylinder 
with two hemispherical endcaps of radius R, we assume that RNase is uniformly distributed on the membrane, and 
we introduce the angle 9 at the hemisphere center between a point on the hemisphere and the long cell axis, see 
Fig. S9. Given an angle element d6, we consider corresponding length element dx along the x axis. We consider 
the amount of RNase along the 9 angle, which is proportional to Rd9, and we project it on the length element 
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FIG. S9: Modeled distribution of RNase on the cell membrane, and its projection on the long cell axis. RNase is uniformly 
distributed on the membrane, and it is represented as an orange layer. The angle element d9 describes an arc on the cell 
membrane (gray) which is then projected on the x axis, where it corresponds to the length element dx (blue). 



FIG. SIO: Steady-state mRNA distributions for the model with 30S and 50S ribosomal subunits and spatially dependent 
mRNA decay, where we fixed the maximal number of allowed translating (T) 30S-50S (70S) pairs per mRNA at n < Umax = 24, 
and we show mRNA species with I < 4 and m < 4 transiently bound (B) 30S and 50S subunits, respectively. (A) mRNA 
profiles for mRNAs with no 30S subunit bound to the translation-initiation site. Total mRNA density ptot{x) (red) and density 
pi,m,nix) of mRNAs with I B 30S subunits, m B 50S subunits, and n 70S pairs (gray),where we show only the curves with I, 
m, and n even for greater clarity. The density po,o,o{x) of ribosome-free mRNAs (green) and the density P 4 , 4 , 24 (a;) of mRNAs 
with the largest number of B, T subunits considered (blue) are also shown. The profiles pi,m,nix), po,o,oix), and ^ 4 . 4 , 24 ( 2 :) are 
normalized to unit area. (B) Same as panel A for mRNAs with a 30S subunit bound to the initiation site. 


dx. As a result, in the hemispherical region £ — R<x<£ we obtain an effective, one-dimensional RNase density 
proportional to where crix) = (2 — is the ratio between the cell cross section 

in the hemispherical endcap and the cross section at midcell. In the central region 0 < x < £ — R, the effective RNase 
density is constant, reflecting the uniform RNase distribution in the surrounding cylindrical cell membrane. The axial 
RNase distribution above will result in an effective, space-dependent mRNA degradation rate /?(a;), which we will 
assume to be proportional to the RNase profile: 


f /3* 0 < X < £ — R, 

Pix) = < £-r<x<£, (S116) 

where the constant has been chosen so as to obtain an average mRNA decay rate (1/f) dx Pix) equal to the 
uniform rate f3 = 3 x 10“^/s of the model discussed in section II. 

In addition, we observe that the hemispherical shape of the cell endcaps results in a reduced cross section of the cell 
perpendicular to the long cell axis, and thus in a smaller available volume for both ribosomes and mRNAs at the cell 
poles. It follows that in the polar region the available volume will be reduced by an amount equal to cix), and we set 
v^ix) v^ix)aix) and ni+m+„(x) (a;)cr(a:). 
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FIG. Sll: Steady-state ribosomal-subunit concentrations and fluxes for the model with SOS and SOS ribosomal subunits and 
spatially dependent mRNA decay. Top: Concentrations C 7 os(a;), C 3 osB(a:), C 3 osf( 2 ;), C 3 os(a;) of 70S, transiently bound (B) SOS, 
free (F) SOS subunits, and of SOS subunits bound to the translation-initiation site, respectively. We also show the concentrations 
Cbosb(x) and C5osf(x) of B and F SOS subunits, respectively. Bottom: Fluxes Jjosix), J3osb(x), J 3 osf(x) of 70S, B and F SOS 
subunits, and flux 7303 ( 2 :) of SOS subunits bound to the translation-initiation site. The fluxes J 5 osb( 2 ;), J 5 osf(x) of B and F 
SOS subunits are also shown. Fluxes are represented along the cell’s long axis depicted in the top panel, the arrow length Is 
proportional to local ribosome flux, and the arrows in the legends correspond to a flux of SO/s. 


Finally, the hemispherical shape of the cell endcaps modifies the binding rates. Indeed, consider a binding rate 
kon and the binder concentration c(x) in our one-dimensional model. First, fcon and c(a:) can be related to the 
three-dimensional binding rate and concentration k^^j^ and C 3 D (x) by equating the number of bindings per unit time: 
konc(x) = k 3 DC 3 B(x). Second, the one- and three-dimensional concentrations are related by c(x) = C 3 D(a:) 7 r(i?sin 0 )^ = 
C3r,(x)a(x)7rR^. Putting everything together, we obtain that the binding rates in the model with hemispherical cell 
endcaps are given by kon = ^on/'^(^)> where = k 3 n,/{'KB?) is the binding rate for a cylindrical cell used in section 
SIOB. 

The mRNA densities, ribosomal concentrations and fluxes, and polysome distributions for the model above with 
spatially dependent mRNA decay are shown in Figs. SIO, Sll, and S12, respectively, where the maximal number of 
allowed 70S pairs per mRNA has been fixed at Umax = 24. 


S12 Results for the late phase of the cell division cycle 

In this section we extend the results from the two-subunit model discussed in section SIO to a cell in the late phase 
of its division cycle. To achieve this, we manually selected 21 dividing cells with a markedly double-lobed nucleoid, 
rescaled their DNA fluorescence profiles to their median cell length of 2^ ~ 4.35/rm, and estimated the nucleoid 
profile along the long cell axis by averaging over multiple cells, see Fig. S13. We chose a density of DNA length 
'^{x) oc 1/{1 -I- exp[^(a:/£ — 2/3)]}{exp[—3C/10(l/2 — x/C)^\ + exp[—3C/10(l/2 -|- x/tj"^]} to reproduce the average 
DNA fluorescence profile in Fig. S13, and we adjusted the model parameters as follows in order to describe the late 
phase of the cell cycle. The parameters that are intensive with respect to the cell length, such as the binding and 
unbinding rates, the diffusion coefficients, and the mRNA decay rate, are taken to be the same as in the model with 
2£ = 3/rm, see section SIOB. On the other hand, extensive parameters were obtained by considering the corresponding 
parameter values of the model with = 3^m, and scaling them linearly with the cell length. As a result, we have 
Aaostot = Wsostot = 8.7 X 10"^, A^mRNA = 7.2 X 10^, WosF = A^sosf = 2% Wot, and Wos = 80% Wot- The nucleoid 
profile a{x) is proportional to <p(x), and its normalization is chosen to achieve the total number of mRNAs above, i.e. 
«tot//3 = iVniRNA- Assuming that both the DNA plectonemic length L and the nucleoid volume V scale linearly with 
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FIG. S12: Polysome distributions for the model with SOS and SOS ribosomal subunits and spatially dependent mRNA decay, 
where we fixed the maximal number of allowed translating (T) 30S-50S (70S) pairs per mRNA at n < rimax = 24, and we show 
mRNA species with I < 4 and m < 4 transiently bound (B) SOS and SOS subunits, respectively. (A) distribution of mRNA 
species with no SOS subunit bound to the translation-initiation site, shown as a heat map of the number 77i,m,n = /q dxpi^m,n{x) 
of mRNAs with I B SOS subunits, m B SOS subunits, and n T 70S pairs. The diameter of each dot is proportional to 77i,m,n, 
where the largest dot in the plot corresponds to Ni^rn,n. ~ 1. (B) Same as panel A for the number dxp'i „^^^{x) of 

mRNAs with a SOS subunit bound to the initiation site, where the largest dot in the plot corresponds to „ ~ 20. 



FIG. SIS: DNA fluorescence along the long cell axis for E. coli cells in the late phase of the division cycle in glucose minimal 
media. Fluorescence for a few representative cells (gray), and resulting average fluorescence over 21 cells in the same phase 
of the division cycle with standard error of the mean (red), both normalized to unit area. Inset: Representative cell near late 
division phase, with ribosomal protein (green) and nucleoid (red). Scale bar: 1 pm. 


the cell length, the normalization of the DNA density (l/i) dxip{x) = LjV is the same as for 21 = 3/rm. 

With this choice of parameters, we solved the reaction-diffusion equations for the two-subunit model by fixing the 
maximal number of allowed 70S pairs per mRNA at rimax = 28: the resulting mRNA densities, ribosomal concentrations 
and fluxes, and polysome distributions are shown in Figs. S14, S15, and S16, respectively. 


S13 Results for different growth rates 

In this section we extend the results from the two-subunit model discussed in section SIO to cells with different 
growth rates. In particular, we consider cells in glycerol minimal and defined rich media, with growth rates ^ 0.56/h 
and 2.3/h, respectively, whose DNA fluorescence profiles are shown in Fig. S17. We will thus consider a density of 
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FIG. S14: Steady-state mRNA distributions in the model with SOS and SOS ribosomal subunits for a long cell in the late 
phase of its division cycle, where we fixed the maximal number of allowed translating (T) 30S-50S (70S) pairs per mRNA at 
n < Wmax = 28, and we show mRNA species with I < 4 and m < 4 transiently bound (B) SOS and SOS subunits, respectively. 
(A) mRNA profiles for mRNAs with no SOS subunit bound to the translation-initiation site. Total mRNA density ptot{x) (red) 
and density pi^m,n{x) of mRNAs with I B SOS subunits, m B SOS subunits, and n 70S pairs (gray), where we show only the 
curves with I, m, and n even for greater clarity. The density po,o,o(x) of ribosome-free mRNAs (green) and the density p 4 , 4 , 28 (x) 
of mRNAs with the largest number of B, T subunits considered (blue) are also shown. The prohles pi^m,n(x), pop,o(x), and 
P4,4,28(x) are normalized to unit area. (B) Same as panel A for mRNAs with a SOS subunit bound to the initiation site. 


DNA length (p(x) to reproduce the profiles in Fig. S17, and adjust all other parameters to the growth conditions along 
the lines of section S12. 

In glycerol minimal media, we selected cells in the mid-phase of the cell division cycle as follows. We considered a cell 
length 2 £ given by ^/^median glycerol ~ ^giucose/^median giucosej where 2 £gi,^cose ~ 3/im is the typical, medium length of 
a cell in the mid-phase of the division cycle in glucose minimal media, see section II. We obtained 2£ ^ 3 p,m, selected 
cells with length within 5% of 2 £, rescaled their DNA fluorescence profiles to a cell length of 2 £, and estimated the 
nucleoid profile along the long cell axis by averaging over multiple cells, see Fig. S17A. Given that £ is close to -^giucosej 
and that Figs. S17A and 2 show that the nucleoid profile in glycerol minimal media is close to that in glucose minimal 
media, all model parameters in glycerol minimal media will be close to those in glucose minimal media, and so will be 
the model predictions. 

Proceeding along the lines of the analysis in glycerol minimal media, for defined rich media we obtained 2 £ ~ 3.5 p,m, 
and we chose a density of DNA length ip(x) oc 1/{1 -|-exp[3/4 ^(a;/£ — 2/3)]}{exp[—C/5(l/2 — x/£Y]-\- exp[—^/5(l/2 -|- 
x/£Y\] to reproduce the observed DNA fluorescence profile, compare Fig. S17B. The intensive model parameters 
are the same as for 2 £ = 2 £giucose = 3 /im, while we estimated the extensive parameters by scaling them with 
respect to the cell length as in section S12. We obtained A^aostot = Nbostot = 7 x 10'^, A^mRNA = 5.8 x 10^, 
N 30 S F = iVsos F = 2 % Wot, Wos = 80 % Wot- In addition, the normalization of the mRNA-synthesis profile a{x) oc '^{x) 
is given by atot = /3WnFi.NA, while the normalization of the DNA density (l/£) fgdx(p(x) = L/V is the same as for 
2£ = 2£giucose- With this choice of parameters, we solved the reaction-diffusion equations for the two-subunit model 
by fixing the maximal number of allowed 70S pairs per mRNA at rimax = 28: the resulting mRNA densities, ribosomal 
concentrations and fluxes, and polysome distributions are shown in Figs. S18, S19, and S20, respectively. 


S14 Model with no transiently bound ribosomes 


In this section we disallow both B 30S and B 50S subunits as follows: First, we set = 0- 

Second, we solve Eqs. (S73) and (S74) by setting Wos = 98% Wot, to take account of the increased number of 70S 
pairs due to the absence of B subunits. 

With this choice of parameters, we solve the reaction-diffusion equations (S90), (S91), (SlOO), and (SlOl). Given 
that there are no B 30S nor B 50S subunits, the only nonzero mRNA densities and mRNA numbers are those with 
I = m = 0, i.e. pofi^n{x), Po,o.n(^), W.o.n and A^o,o,n- Fig- ^21 we show the resulting mRNA profiles, the total 
mRNA densities, and the mRNA numbers, for both mRNAs with and without a 30S subunit bound to the translation- 
initiation site, while in Fig. S22 we show the concentrations and fluxes of ribosomal subunits. Overall, these results 
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FIG. S15: Steady-state ribosomal-subunit concentrations and fluxes in the model including SOS and SOS ribosomal subunits 
for a long cell, in the late phase of its division cycle. Top: Concentrations 0703 ( 2 :), C 3 osB(a:), C 3 osF(a:), € 303 ( 2 ;) of 70S, transiently 
bound (B) SOS, free (F) SOS subunits, and of SOS subunits bound to the translation-initiation site, respectively. We also show 
the concentrations 05033 ( 2 ;) and 05033 ( 2 ;) of B and F SOS subunits, respectively. Bottom: Fluxes J 7 os( 2 ;), J 303 b( 2 ;), J 3 osf( 2 ;) of 
70S, B and F SOS subunits, and flux Jsosi^) of SOS subunits bound to the initiation site. The fluxes J 503 b( 2 :), J 5 osf( 2 ;) of B 
and F SOS subunits are also shown. Fluxes are represented along the cell’s long axis depicted in the top panel, the arrow length 
is proportional to local ribosome flux, and the arrows in the legends correspond to a flux of SO/s. 



FIG. S16: Polysome distributions in the model including SOS and SOS ribosomal subunits for a long cell in the late phase 
of its division cycle, where we fixed the maximal number of allowed translating (T) SOS-SOS (70S) pairs per mRNA at n < 
n-max = 28, and we show mRNA species with I < 4 and m < 4 transiently bound (B) SOS and SOS subunits, respectively. (A) 
distribution of mRNA species with no SOS subunit bound to the translation-initiation site, shown as a heat map of the number 
Ni^rn.n = dxpi^rn,n(x) of mRNAs with I B SOS subunits, m B SOS subunits, and n T 70S pairs. The diameter of each dot 
is proportional to Ni^rn.n, where the largest dot in the plot corresponds to W,m,n ~ 2. (B) Same as panel A for the number 
~ fa dxpi „^ „(x) of mRNAs with a SOS subunit bound to the initiation site, where the largest dot in the plot corresponds 
to ~ 30. 
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FIG. S17: DNA fluorescence along the long cell axis for E. coli in glycerol minimal and defined rich media. (A) Glycerol 
minimal media: fluorescence for a few representative cells (gray) and resulting average fluorescence over 81 cells with standard 
error of the mean (red), both symmetrized and normalized to unit area. Inset: representative cell grown in glycerol minimal 
media, with ribosomal protein (green) and nucleoid (red). Scale bar: 1 /rm. (B) Same as panel A for defined rich media, averaged 
over 43 cells. 



FIG. S18: Steady-state mRNA distributions for the model including SOS and SOS ribosomal subunits for cells in defined rich 
media with growth rate ~ 2/h, compare Fig. S17B, where we fixed the maximal number of allowed translating (T) 30S-50S 
(70S) pairs per mRNA at n < rimax = 28, and we show mRNA species with I < 4 and m < 4 transiently bound (B) SOS and 
SOS subunits, respectively. (A) mRNA profiles for mRNAs with no SOS subunit bound to the translation-initiation site. Total 
mRNA density ptot(x) (red) and density pi,m,n{x) of mRNAs with I B SOS subunits, m B SOS subunits, and n 70S pairs (gray), 
where we show only the curves with I, m, and n even for greater clarity. The density po,o,o(x) of ribosome-free mRNAs (green) 
and the density P 4 , 4 , 28 (®) of mRNAs with the largest number of B, T subunits considered (blue) are also shown. The profiles 
pi,m,n(x), po, 0 , 0 ( 2 ;), and P4,4,28(2;) are normalized to unit area. (B) Same as panel A for mRNAs with a SOS subunit bound to 
the initiation site. 


confirm those including B subunits. An important feature distinguishing the existence from the absence of B ribosomes 
is the ribosome efficiency: without B ribosomes, this efficiency is e = 2.39 x I0“^/s, which is ~ 23% higher than the 
efficiency including B ribosomes, compare Section SIOD. Put simply, as the total number of ribosomal subunits is 
conserved in this comparison, the inclusion of transient ribosomal binding reduces the number of 30S and 50S subunits 
that can form translating 70S pairs, thus lowering the ribosome efficiency. In addition. Fig. S21 shows that there is 
an average number of ^ 12 subunits, i.e. 30S and 70S, per mRNA. This average mRNA occupancy is smaller than 
that obtained including B subunits in section SIOD: The absence of B ribosomes thus implies that excluded-volume 
effects are slightly weaker, and that the dip in the total mRNA profiles in Fig. S18 is no longer present in Fig. S21. 
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FIG. S19: Steady-state ribosomal-subunit concentrations and fluxes in the model including SOS and 50S ribosomal subunits 
for cells in defined rich media with growth rate ~ 2/hr, compare Fig. S17B. Top: Concentrations Cjos{x), C 3 osb( 2 ;), C 3 osF(a:), 
C 3 os(*) of fOS, transiently bound (B) SOS, free (F) SOS subunits, and of SOS subunits bound to the translation-initiation site, 
respectively. We also show the concentrations C 5 osB(a:) and C 5 osf( 2 ;) of B and F SOS subunits, respectively. Bottom: Fluxes 
^ 703 ( 2 :), >73osb(2:), J 3 osf( 2 ;) of 70S, B and F SOS subunits, and flux of SOS subunits bound to the initiation site. The 

fluxes J 5 osb{x), J 5 osF(a;) of B and F SOS subunits are also shown. Fluxes are represented along the cell’s long axis depicted 
in the top panel, the arrow length is proportional to local ribosome flux, and the arrows in the legends correspond to a flux of 
SO/s. 



FIG. S20: Polysome distributions for the model with SOS and SOS ribosomal subunits for cells in defined rich media with 
growth rate ~ 2/h, compare Fig. S17B, where we fixed the maximal number of allowed translating (T) S0S-50S (70S) pairs per 
mRNA at n < Umax = 28, and we show mRNA species with Z < 4 and m < 4 transiently bound (B) SOS and SOS subunits, 
respectively. (A) distribution of mRNA species with no SOS subunit bound to the translation-initiation site, shown as a heat 
map of the number Ni^rn,n = Jq dxpi^rn,n{x) of mRNAs with I B SOS subunits, m B SOS subunits, and n T 70S pairs. The 
diameter of each dot is proportional to Ni^rn,n, where the largest dot in the plot corresponds to W,m,n ~ 1. (B) Same as panel 
A for the number ,j(a:) of mRNAs with a SOS subunit bound to the initiation site, where the largest dot in 

the plot corresponds to ~ 20. 
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FIG. S21: Steady-state mRNA distributions for the model including SOS and SOS ribosomal subunits, with no transiently bound 
(B) SOS and SOS subunits. The maximal allowed number of 70S pairs per mRNA is rimax = 28. (A) mRNA profiles for mRNAs 
with no SOS subunit bound to the translation-initiation site. Total mRNA density ptot{x) (red) and density pop,nix) of mRNAs 
with n 70S pairs (gray). The density po,o,o(x) of ribosome-free mRNAs (green) and the density po,o, 28 (x) of mRNAs with the 
largest number of 70S pairs considered (blue) are also shown. The profiles po^o,n(x), Po,o,o(x), and po,o, 28 (®) are normalized to 
unit area. Inset: distribution No,o,n = fg dxpo,o,n(x) of mRNAs with n 70S pairs in the right half of the cell as a function of n. 
(B) Same as panel A for mRNAs with a SOS subunit bound to the initiation site. 
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FIG. S22: Steady-state ribosomal-subunit concentrations and fluxes for the model including SOS and 50S ribosomal subunits, 
with no transiently bound (B) SOS and SOS subunits. Top: Concentrations 0703 ( 2 ;), C 3 osF(a;) of 70S and free (F) SOS subunits, 
concentration 0303 ( 2 :) of SOS subunits bound to the translation-initiation site, and concentration C 5 osf( 2 :) of F SOS subunits, for 
the right half of the cell. Bottom: Fluxes J 7 os( 2 ;), J 3 osf( 2 :) of 70S and F SOS subunits, flux Jsosix) of SOS subunits bound to 
the initiation site, and flux Jsospix) of F SOS subunits. Fluxes are represented along the cell’s long axis depicted in the top 
panel, the arrow length is proportional to local ribosome flux, and the arrows in the legends correspond to a flux of S/s. 


S15 Force balance between nucleoid and polysomes determines nucleoid size 

In this section we will calculate the force exerted by polysomes and F ribosomes on the nucleoid as well as the 
entropy of compaction of the nucleoid, and the resulting force-balance equation will allow us to determine the nucleoid 
size. 

For the sake of clarity, let us first focus on F ribosomes. We compute the total force exerted by the nucleoid on F 
ribosomes, and the total force by F ribosomes on the nucleoid will necessarily be equal and opposite. To begin, we 
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observe that the F-ribosome current (S3) is the sum of two terms which can be interpreted as follows: The first term, 

(S117) 

ax 

is the current due to the Brownian diffusion of F ribosomes from high- to low-concentration regions, and it is thus 
proportional to the concentration gradient, with a space-dependent diffusion coefficient D'pv^{x). The second term, 

Junhix) = D-pc^{x) —(S118) 

is the current due to the force locally exerted by the nucleoid on F ribosomes; it is proportional to the gradient of 
Uf(x), thus it represents the tendency of the nucleoid to push ribosomes from regions with lower available volume 
to regions with higher available volume. The current JDNA(a^) is thus the result of a net force f{x) applied by the 
nucleoid on each F ribosome at position x, with the force and the current related by 


JDNA(a;) = n{x)f{x)cj.{x), 

where the mobility /x(x) is related to the Brownian diffusion coefficient DpVp{x) by the Einstein relation 


(S119) 


DpVf{x) = n{x) ke T, 


(S120) 


where ke is Boltzmann’s constant and T is the temperature [25]. From Eqs. (S118), (S119), and (S120) we obtain 
~ ■ Multiplying f{x) by the local F ribosome concentration, integrating with respect to x, and 

reversing the sign of the force, we obtain the total force exerted by F ribosomes in the right half of the cell on the 
nucleoid, which reads —ke T dx ■ Finally, the combined force exerted by both F ribosomes and polysomes 

on the nucleoid is obtained by incorporating the mRNA concentrations Pm,n{x) and available volumes Vm+n(x) in the 
expression above, and summing over all mRNA species. We obtain 




= —ke T f dx 
Jo 

= —ke T f dx 
Jo 
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(S121) 


konCFjx) 
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1 dvp{x) 
v-p{x) dx 


Pn{x) 


where the subscript ‘in’ denotes a force directed towards the inner part of the cell, and in the second line we used Eqs. 
(S20) and (S25), which are valid in the rapid-equilibrium limit, to compute the sum with respect to m. Note that J^in 
depends on the DNA profile (p{x) through the available volumes Uf]^) and Vr{x), compare Eqs. (S9) and (Sll). 

The force exerted by F ribosomes and polysomes results in the compaction of the nucleoid [8, 26]. At mechanical 
equilibrium, this force must balance the “spring” force exerted by the compressed nucleoid, which can be estimated 
in terms of the entropy of a self-avoiding DNA polymer. We approximate this entropy S by treating the plectonemic 
chromosomal DNA as a sequence of N joined segments, each with length confined in a volume V, yielding the 
confinement dependence of the entropy ]27] 


^ B 2 ylj-i > 


(S122) 


where ^ = 200 nm is twice the estimated persistence length of the plectoneme [8], the total number of monomers 
N = Lj^ = 7.5 X 10^ is obtained as the ratio between the total plectoneme length L and the segment length and 


= 0.592 (S123) 

is the exponent describing the mean end-to-end distance for a self-avoiding polymer ]27]. Given that the exponent 
(S123) results from a numerical simulation where the polymer is represented as a self-avoiding walk on a cubic lattice 
[27], the estimate (S122) applies to self-avoiding polymers whose thickness is comparable to the segment length. It 
follows that for a DNA plectoneme the expression (S122) should be, in principle, modified to account for the difference 
between the thickness and segment length. However, in what follows we will show that Eq. (S122) quantitatively 
reproduces the observed nucleoid size, thus showing that the thickness correction above is negligible. 

Proceeding along the lines of Eq. (S7), we assume that the DNA is confined in a cylindrical region at midcell: We 
denote the length of this region by 2xo, and the DNA profile is ip(x) <x. l/{l-|-exp]C(a; —xq)/^]}, where p is normalized 
according to the condition (S8) on the total DNA length. The volume V in Eq. (S122) can thus be estimated as 
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FIG. S23: Steady-state mRNA and polysome distributions for the model including the force balance between nucleoid and 
polysomes. Total mRNA density ptot{x) (red) and density pm,n{x) of mRNAs with m transiently bound (B) ribosomes and n 
translating (T) ribosomes for 0 < m < mmax = 8 and 0 < n < rimax = 24 (gray). The density pofi(x) of ribosome-free mRNAs 
(green) and the density pg, 24 (a:) of mRNAs with the largest number of T and B ribosomes considered (blue) are also shown. 
The profiles pm,n(x), pofi(x), and pg, 24 (a:) are normalized to unit area. Inset: distribution of mRNA species, shown as a heat 
map of the number of mRNAs with m B ribosomes and n T ribosomes in the right half of the cell. The maximal number 

of T ribosomes per mRNA used in our model, rimax = 24, is indicated. 


V = 2'kR^xq, and the outward entropic force exerted by the nucleoid on F ribosomes and polysomes in each half of 
the cell is obtained from the derivative of the entropy as 




2 Oxq 
2 

liB T -=-xr , 

6 (27ri?2)2/3a.5/3 


(S124) 


where the factor 1/2 in the first line results from the fact that we compute the outward force in one half of the cell. 
In what follows we will thus use Eqs. (S121) and (S124), combine the resulting force-balance equation 




Flmt = 0 


(S125) 


with the reaction-diffusion Eqs. (S121) and (S124), and solve self-consistently for Pn{x), CF(a:), and for the nucleoid 
size xq. The results are shown in Figs. S23 and S24, where we show mRNA concentrations and polysome distributions, 
the nucleoid profile p{x), the calculated value of « 0.5, as well as the F ribosome concentrations and fluxes. 

Finally, we provide a simple analytical estimate for the nucleoid size as a function of the total number of mRNAs, 
and other relevant parameters. Given that the ribosome-to-niRNA ratio is much larger than unity, compare section 
II, we assume that mRNAs at the poles are fully excluded from the nucleoid by virtue of having a sufficient number 
of bound ribosomes: we thus approximate the nucleoid profile in one half of the cell as a step function with the step 
&t X = xq, i.e. Vm+n(x) is equal to zero for 0 < a; < xo and to unity for xq < x < 1. Given that the size of an F 
ribosome is significantly smaller than that of a typical polysome, F ribosomes can easily penetrate the nucleoid, thus 
the pressure that they exert on the nucleoid is much smaller than the force that polysomes exert on the nucleoid: 
Indeed, in Eq. (S121) we denote the contribution to F/n given by F ribosomes and mRNAs by F/nP s-nd F/nmRNA> 
respectively, and for the reference conditions of Figs. 2-4, we have -FinF/i-^inF +F'inmRNA) 2 x 10“^. We thus 
neglect the contribution by F ribosomes in the first line of Eq. (S121), and we obtain 


^^max ^max 




771—0 77 — 0 


n dx 




(S126) 


where in the first line we used the fact that at equilibrium pm,n{x)/Vm+n{x) is independent of x, compare section S2, 
and we denoted by the concentration of mRNA species m, n outside the nucleoid, and by 


out ^ -^mRNA 

2{^-Xo) 


(S127) 
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FIG. S24: Nucleoid profile, steady-state ribosome concentrations, and ribosome fluxes for the model including the force balance 
between nucleoid and polysomes. Top: DNA density ‘fi(x) normalized to unit area and concentrations Ct(x), Cb(x), and Cf(x) 
of translating (T), transiently bound (B), and free (F) ribosomes in the right half of the cell (compare Fig. 1). The nucleoid 
size xo/£ is also marked. Bottom: Fluxes of T, B, and F ribosomes along the cell’s long axis depicted in the top panel. The 
arrow length is proportional to local ribosome flux, and the arrows in the legends correspond to a flux of 20/s. 
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FIG. S25: Nucleoid size as a function of number of mRNAs, as predicted by force balance between nucleoid and polysomes. We 
plot the ratio xo/£ between the nucleoid width and cell length as a function of the total number of mRNAs in the cell, from the 
solution of Eq. (3), which relies on the simplifying assumption that mRNAs cannot penetrate the nucleoid. The mRNA number 
NmRNA = 5 X 10^ used for our analysis in glucose minimal media and the predicted nucleoid size xoft ~ 0.5 are indicated with 
a dot. 


the total mRNA concentration outside the nucleoid. Combining Eqs. (S124), (S125), (S126), and (S127) we obtain 
Eq. (3), which can be solved numerically for the nucleoid size xq. The resulting values for the nucleoid size Xq are 
shown in Fig. S25, where we plot xo/£ as a function of the total mRNA number Ai„iRNA keeping all other parameters 
in Eq. (3) fixed. 
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